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ABSTRACT 

We consider the high radio frequency (15 GHz - 353 GHz) properties and variability 
of 35 Brightest Cluster Galaxies (BCGs). These are the most core-dominated sources 
drawn from a parent sample of more than 700 X-ray selected clusters, thus allowing 
us to relate our results to the general population. We find that >6.0% of our par¬ 
ent sample (>15.1% if only cool-core clusters are considered) contain a radio-source 
at 150 GHz of at least 3mJy (wlxlO 23 W Hz~ x at our median redshift of z«0.13). 
Furthermore, >3.4% of the BCGs in our parent sample contain a peaked component 
(Gigahertz Peaked Spectrum, GPS) in their spectra that peaks above 2 GHz, increas¬ 
ing to >8.5% if only cool-core clusters are considered. We see little evidence for strong 
variability at 15 GHz on short (week-month) timescales although we see variations 
greater than 20% at 150 GHz over 6-month timesframes for 4 of the 23 sources with 
multi-epoch observations. Much more prevalent is long-term (year-decade timescale) 
variability, with average annual amplitude variations greater than 1% at 15 GHz be¬ 
ing commonplace. There is a weak trend towards higher variability as the peak of the 
GPS-like component occurs at higher frequency. We demonstrate the complexity that 
is seen in the radio spectra of BCGs and discuss the potentially significant implications 
of these high-peaking components for Sunyaev-Zel‘dovich cluster searches. 

Key words: radio continuum: galaxies - clusters: general 


1 INTRODUCTION 

Energetic feedback from accreting supermassive black holes 
* E-mail: m4hogan@uwaterloo.ca (MTH) (SMBH) is now widely accepted to play an integral role in 
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the formation and evolution of Universal structure. Such 
action is commonly invoked to explain a variety of phenom¬ 
ena such as the hi gh-end curtailment of the galaxy lumi¬ 
nosity fu nctio n (e.g. Benson et a l. 20031 ; [Bower et al. ] 20061: 


Croton_e£_atJ_ 20061'). t h e M bh vs M; , u iae correlation (e.g 

Magorrian et al. 11 19981 : ISilk fc Rees I Il998l l and the sym¬ 


metry seen in the cosmic histories of both star formation 
and active galactic nuclei (AGN) (e .g. iMerloni et al. I [2004 
ISpringel et al. Il2005l : I T I o t> k i n si 1 2 012T) . 

In the most massive systems, at the centres of galaxy 
clusters, it is now well established that the AGN-action 
of the centrally located Brightest Cluster Galaxy (BCG) 
prevent s runaway coo ling (the classical Cooling Flow Prob¬ 
lem : see IFabian II 1994 . for a review) by imparting energy to 
its surroundings thr ough mechan ical feedback (see the re- 
cent reviews by iMcNamara fc Nulsen 1 120071 . 120121 : iFabian I 
l2012fl . Such a ctions are invoked to explain the defi cit of 
cold gas (e.g. lEdge I l200ll: ISalo me fc Combes 1 12003 ) and 


star formation ( O’Dea et al. I ~ 20081 : Rafferty et al. 


200811 


compared to what would be expected from cooling dom¬ 
inated systems, as well as t he dearth of gas at interme- 
diate cooling temperatures dPeterson et al. I l2003r) . Strik¬ 
ing observational evidence for this mechanical feedback 
is seen through the presence of X-ray cavities in the 
intra-cluster medium (ICM) (e.g. McNmnaraef_at 200d: 


IFabian et al. l20~od : lHlavacek-Larrondo et al. 112012a! '). These 

cavities, which are inflated by the actions of co-spatially 
observed radio-emitting plasma, subsequently re-distribute 
accretion energy from the central super-massive black hole 
(SMBH) to the surroundings as they buoyantly rise. 

There has apparently been a fine balance between 
heating and cooling within clusters in place for at least 
half the Hubble time (e.g. IVikhlinin~l 120061 : iPratt et al. I 
l20ld : lMcDonaId~ et al. ]»). Furthermore, the energy im¬ 
parted by AGN activity doe s appear to be sufficient to 
offset cooling on avera ge (e.s. IMcNamara fc Nulsen~ll2007l ; 
iDunn fc Fabianl l2008l h However, this is not the case at 
all times suggesting that periods of cooling must be inter¬ 
spersed by periods of AGN energy injection. Added to this 
is the growing realisation that in the most settled cool-core 
(CC) clusters where there is a central peak in the cluster 
X-ray surface profile indicative of substantial cooling, the 
BCGs have a radio-loud duty cycle approaching unit y (e.g. 

I Burns iFl 99()l : iMittal et al. 1120091 ; IHogan et al. Il2015al '). This 
suggests that in these systems there must be cyclic activ¬ 
ity whereby the BCG is more active in some periods than 
others. 

Radio observations are integral for tracing this mechan¬ 
ical feedback. Ho wever, most studie s have been carried out 
at <1.4 GHz fe.g. iBest et al. 1120071 1 and so the radio prop¬ 
erties of BCGs in the ~10-300 GHz range are somewhat 
poorly constrained. Single-dish surveys of the sky at the 
higher end of this ra nge are typically shallow with relatively 
low resolution (e.g. IPlanck Collaboration 1 |2014a) whereas 
even at the lower end, small beam sizes make interferomet¬ 
ric surveys of any sizeable area both difficult and expensive. 
However, several recent surveys have allowed for huge ad- 


10 GHz fe.g. 10C at 15.7 GHz, Franzen et al. 201ll: AT20G 

at 20 GHz. iMurohy et al. 201C 

and the AT20-deep (pilot) 

also at 20 GHz. Massardi et al. 

201 ll: Franzen et al. 2014|j. 


Consequently, only a few of the brightest BCGs have well 
charact erised radio-spectr a in this crucial spectral range. 

In IHogan et al. I d2015al l (hereforth H15a) we consid¬ 
ered the radio properties of BCGs in a parent sample of 
over 700 X-ray selected clusters comprising the BCS, eBCS 
and REFLEX cluster catalogues llEbeling et al. Ifl998l . l2000l : 
iBohringer et al 20041. respectively). These catalogues are 
X-ray flux limited, hence our sources are not selected on 
radio priors and should be representative of the general 
cluster population. The clusters in the parent sample were 
split into CCs and non cool-cores (NCCs) using the pres¬ 
ence of optical emission lines around the BCG. Such lines 
are only found in systems with central cooling times less 
than 5xl0 8 years, equivalent to a central entropy less than 
30 keV cm 2 and hence can be use d as a proxy for the 
cluster state llC ayagnol o et al. 1 120081 : iRaffertv et al. 1 120081 : 
ISanderson et al 20091 '). The radio spectral energy distribu¬ 
tions (SEDs) of the BCG in 246 of these were populated, 
typically between 74 MHz and 20 GHz, and decomposed 
into active and inactive components attributable to ongoing 
and historical accretion respectively (see H15a). Not only 
was the radio duty cycle of BCGs in CCs seen to be substan¬ 
tially higher than in NCCs, it was found that the majority of 
CC-hosted BCGs showed evidence for ongoing core activity 
that manifests itself primarily as a spectral flattening above 
a few GHz. Often this emission is missed in low frequency 
surveys. Further confounding the lack of information, in¬ 
creased variability is postulated to higher radio frequencies 
as the emission is expected to originate from increasingly 
smaller physical scales. 

In this paper we select a sub-sample of the BCGs stud¬ 
ied in H15a, believed to contain the most active cores. These 
sources all reside in CC clusters where active feedback is 
prevalent. Furthermore, the pre-eminence of the radio core 
component in the radio SEDs of these sources indicates that 
the SMBH is actively accreting at a significant rate. We have 
observed them with a variety of facilities to extend their ra¬ 
dio coverage up to 353 GHz as well as observing a number of 
them at different epochs which, alongside historical observa¬ 
tions, allows us to study their variability. We re-iterate that 
these sources are from an X-ray selected cluster catalogue. 
By selecting the sources with the highest expected core flux 
we are able to search for variability via short observations. 
This permits us to explore the origin of the point-like central 
radio emission as well as to constrain the amplitude of vari¬ 
ation in the accretion rate during active periods of ongoing 
feedback. Whilst the feedback powers derived from X-ray 
cavities trace the AGN energy output averaged over tens 
of megayears, this shorter term variability provides insights 
into the more instantanous processes within the core. 

One system that has been well-monitored at radio fre¬ 
quencies greater than a few GHz, and indeed constitutes 
one of the most well-studied examples of active AGN feed¬ 
back in a cl uster core, is NGC1275/3C84 in the Perseus 
Cluster (e.g. Bohringer et al. 19931 : IConselice~ et al. 200ll : 
lAbdo et al. 1120091 4. Large amplitude variations in the radio 
spectrum of this source have been known for many years 


I Paulinv-Toth & KeUermannJ 


196' 


Recently Dutson et al. (I2014T) undertook a compre¬ 


hensive study of the radio and gamma-ray properties of 
NGC1275, considering its radio variability over five decades 
in both time and frequency. The radio spectrum consists 
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of a steep spectrum power-law at frequencies below ap¬ 
proximately 1 GHz and an inversion above this leading 
to a peaked profile. I t should be noted that the source is 
not s trongly beamed jKrichbaum et fflZ~lll992l : iNagai et al. I 
l2Q10h . The power-law component, attributed largely to 
t he presence of extended lobes and a 300kpc mini-halo 
dBurns et al. I Il992h . is constant in its flux. However the 
peaked component is found to vary significantly in both flux 
(more than an order of magnitude) and turnover frequency 
on few year/decade timescales. Such variations have been 
previously linked to individual components in the jet on mil- 
liarcsecond scales, as recovera ble using Very Long Baseline 
Interf erometry (VLBI, e.g . ISuzuki et al. I 20121 1. Interest¬ 
ingly, iDutson et al. I d2014l l find compelling correspondence 
between this few-year variation of the high radio-frequency 
peaked component and the high energy gamma-ray emission 
but no strong connection between the short-term ‘flaring’ 
seen in the gamma-rays and the 1.3mm hux. In a study o f 
the core X-ray properties of 57 BCGs, iRussell et al. I (|2013l l 
found that roughly half contained an X-ray point source at 
Chandra resolution. It is worth noting that three of these 
(A2052, Hydra-A and M84) were seen to vary over similar 
6-month to decade timescales, similar to the radio emission 
in NGC1275. 


One of the aims of the current paper is to investigate 
whether the high radio frequency properties of NGC1275 
mark it out as a peculiar object or whether such periods of 
high activity in the spectral range above 10 GHz are common 
amongst the BCG population. 

Re cent ly the Sunyaev-Zel‘dovich (SZ) effect 
llSunvaev &; Zel‘dovich I Il972i f has bee n used to compile 
large c atalogues of galax y clust er s (e.g. VaijHeriindg__eLfflLJ 


20 id: Reichardt_ei_al. I 20131: Marriage et al. I 1201 lal: 


HAsselheld et all 20131 : iPlanck Collaboration I 2011al . 20131 . 
201461 1 . Unresolved radio sources present a significant sys¬ 


tematic for these searches and may lead to underestimated 
or completely re moved SZ de c remen t in the 15-20 0 GHz 
range (see e.g. i Knox et al. I |2004 ICoble et al. I 120071 : 
iLin fe Mohr I l2009l l . Furthermore, the single-dish nature 


of many SZ-observatories means that often the removal 
of contaminating point sources has to rely on uncertain 
extrapolation of higher resolution but lower frequency 
data. Our results therefore have potentially significant 
implications for SZ studies of clusters. 

This paper is arranged as follows. In Sections [2] and [3] 
we describe the sample selection, data collection and reduc¬ 
tion. Our results are presented in Section [4] and discussed in 
Section [5] before we draw conclusion in Section [6] We have 
used a standard ACDM cosmology with: 9 m = 0.3, 9 a = 
0.7, Ho = 70 km s -1 Mpc -1 . We use the spectral index con¬ 
vention S oc v~ a . Unless other stated, we use the name of 
the parent cluster to refer to its BCG. 


2 SAMPLE 

The sample of sources chosen for this study were selected 
primarily from H15a as having the brightest (> lOmJy at 
5 GHz), flat-spectrum cores (a < 0.5) so a detection above 
100 GHz was possible. The H15a sample covers an all-sky, 
X-ray flux-limited sample of over 700 clusters as outlined 
above, but the number of sources matching these flux and in¬ 


dex cuts is relatively small (< 30 or < 4%) . To increase the 
target list we added seven bright (>50mJy at 5 GHz) sources 
either in fainter clusters and/or clusters mis-identihed un¬ 
til n ow. Four of these sources ar e in fainter clusters 
(All. lPerlman et al. Il99£ ; 4C+55.16, I wasawa et al. 19991 : 
A227o7 lHealev et al. Il2007 : and RXJ2341+00) and three in 


clusters above the eBCS/REFLEX flux limit (RXJ1350 +09; 


RXJ1832+68. Bohrineer et al. 

200jQ!; Gioia et al. 20031; and 

E1821+64. Russell et al. 201C 

). All seven of these sources 


have been previously identified as AGN, given the associa¬ 
tion of radio and X-ray emission. However, only E1821+64 
is actually strongly (more than 50%) contaminated in the 
X-ray by an AGN, in this case a QSO. In the other six cases 
the cluster has previously been mis-identihed as a BLLac 
as a result of them having seemingly hat radio spectra and 
so the X-rays had been attributed to the AGN. However, 
all six of these clusters have central galaxies with strong, 
narrow optical line emission, characteristic of cooling How 
BCGs llCrawford et al. Ill999l l and in all other aspects are 
similar to the sources selected in H15a. Therefore, we pro¬ 
pose a re-identification of each source such that the X-rays 
are predominantly from the cluster and not a central AGN. 

The source selection is by no means complete but is rep¬ 
resentative of the brightest, core-dominated radio sources in 
cluster cores. Therefore we believe that the spectral and vari¬ 
ability properties we determine for this radio-bright sample 
can be used to constrain the properties of the complete, X- 
ray selected sample as a whole. 


3 DATA 
3.1 GISMO 

GISMO is a 150 GHz (2mm) bo lometer camera built by 
the Goddard Space Flight Centr e (IStaguhn et al. 20081) for 
use on the IRAM-30m telescopqj, which is located at an 
altitude of 2850m on Pico Veleta, Spain. GISMO operates as 
a visitor instrument, being operable for around a two week 
period every six months. During its time on the telescope 
GISMO is operated in a shared risk pool observing mode, 
during which telescope focus and pointing observations are 
regularly performed on IRAM calibrator sources. 

We obtained data from 3 epochs, using GISMO to ob¬ 
serve 29, 24 and 17 sources in the April 2012, November 
2012 and April 2013 observing runs respectively, with as 
many source overlaps between runs as possible (see table 
H. Sources were typically observed for 5-10 minutes using 
a standard lissajous scanning script. The telescope has a 
FWHM of 16.7” and the absolute hux calibration of GISMO 
is found to have a typical 8% uncertainty. 

Data were reduced using the GISMO spec ihc section 
of the CRUSH software package iKovacs 2008! ). Pointing 
and hux models specific to each run are updated during 
each GISMO run and the package itself is being continually 
updated. We therefore reduced our data using a version of 


1 IRAM is supported by INSU/CNRS (France), MPG (Germany) 
and IGN (Spain) 

2 Performance reports available for each run at 
http: //www.iram.es/IRAMES/mainWiki/Continuum/GISMO/M 
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CRUSH (2.15-2) that post-dates all of our data and hence 
contains optimised parameters for each of our data-sets. 

Most of our sources have signal-to-noise ratios (SNRs) 
less than 10 for each scan. We therefore found the -faint op¬ 
tion within CRUSH to deliver the best results in most cases. 
For non-detections we re-mapped the data using the -deep 
filtering. Only one source (A1885) was recovered using -deep 
but not -faint. Since the -deep option is known to over-filter 
sources with SNR greater than 5 and will create negative 
flux haloes, we took our flux measurements from the maps 
created using the -faint option. A1885 was recovered with 
a SNR of «4 and hence measuring the flux from the -deep 
map for this source is not believed to introduce any addi¬ 
tional error. The specific filtering mode used for each source 
is shown in table [5] Atmospheric extinction is automatically 
corrected for within CRUSH, using an estimate of the most 
recent tau225 GHz value as measured by the on-site tau- 
meter. However, some of our data (particularly during the 
Nov. 2012 run) suffered from a faulty tau-meter not inserting 
the correct measurement into the FITS file. Time-dependent 
records of the on-site tau values were recorded manually dur¬ 
ing each run, which were compared to the values recorded 
automatically in the log files. To ensure consistency, we man¬ 
ually inputted the best estimate zenith tau-value using the 
-tau.225 GHz option for each of our source reductions. 

Maps were produced and fluxes extracted using the 
CRUSH tool ‘show’. Gaussian fits were performed for each 
source with the resultant flux and FWHM reported along 
with source peak and map RMS in tabled Only Hydra-A 
showed significant extent and hence its flux was extracted 
from a user-defined region. Fluxe s were additi onally verified 
using the kvis data-display tool dGooch II 1996T) . 

3.2 CARMA 

Twenty three sources were observed at 9 0 GH z using the 
CARMA interferometer fe.g. lWoodv et a/~ll2004l ~l in D-array 
between 21st May - 15th June 2012, of which twenty over¬ 
lap with our GISMO sources. These observations were per¬ 
formed in queue mode by initiating pre-determined blocks, 
which could be started and abandoned part way through, de¬ 
pending on observing conditions and at the discretion of the 
telescope operator. Each block contained an observation of 
a planet for primary flux calibration. This observation is run 
through CARMA pipelines to provide a primary flux scale 
for the observations. Additionally, a strong bandpass cali¬ 
brator was observed for each block (the bandpass calibrator 
used was either 3C279, 3C345, 3C454.3 or 3C446, depending 
on observability). Science targets were visited several times 
at various hour angles to maximise uv-coverage and sand¬ 
wiched between 1 minute observations of a nearby phase 
calibrator. 

Data were reduced by the standard methods using the 
CARMA optimised version of the MIRIAD data-reduction 
software package (developed and maintained by University 
of Maryland). Some automatic flagging of bad data is per¬ 
formed, but the visibilities were also inspected and flagged 
interactively. Files with basic calibration, taken at the time 
of observation, were provided and used as a calibration 
‘starting point’. Corrections for antenna positions were per¬ 
formed using the most co-temporal antpos models provided 
by the observatory. Further bandpass, phase and flux cor¬ 


rections were applied using standard MIRIAD tasks as de¬ 
termined from the observed calibration sources. 

Deconvolution and cleaning of the subsequent maps 
were also performed using MIRIAD and the kvis data- 
display tool used to extract fluxes, verified by additionally 
extracting fluxes using the AIPS task JMFIT. 


3.3 AMI 

The Arcminute Microkelvin Imager (AMI, IZwart et al. I 
l2008h is a versatile instrument located at the Mullard Ra¬ 
dio Astronomy Observatory (MRAO). Optimised for study 
of the SZ effect, AMI consists of two interferometer arrays: 
the Small Array (SA), with ten 3.7 m antennas with base¬ 
lines of 5-20 m; and the Large Array (LA) with eight 12.8 m 
antennas with baselines of 18-110 m. These arrays are thus 
optimised for large (3 - 10 arcmin) and small (30 arcsec 
- 3 arcmin) scale observations respectively. This set-up al¬ 
lows study of the large-scale SZ effect with accurate char¬ 
acterisation of the contaminating radio source environment. 
Both arrays measure the I+Q polarisation at the central fre¬ 
quency of 16 GHz. Data are taken in six channels, each of 
width 0.72 GHz, over the range 13.9-18.2 GHz, allowing a 
local spectral index to be determined for strong sources. 

We used the AMI-LA to observe seventeen of our 
sources, with each target visited either two or three times 
in 2012 (see Table [3j) . Observations consisted of 8 minute 
integrations sandwiched between 1 minute phase calibrator 
scans. In the April-June scans sources were themselves used 
for self-calibration. For the September follow-ups we addi¬ 
tionally observed nearby phase calibrator sources. 3C 48 and 
3C 286 were observed for use as amplitude calibrators. 

Data reduction was performed using the local in-house 
software Reduce. Reduction pipelines were used to apply 
amplitude and phase calibration, flag for telescope errors 
and Fourier transform to produce frequency channels. Ad¬ 
ditional flagging for bad data was carried out interactively 
through visual inspection of the channel data. 

Data are written out of Reduce as uv-FITS files which 
were read into AIPS for deconvolution and cleaning using 
IMAGR. Maps were produced for each channel in addition 
to a single image for the full bandwidth. None of our sources 
were found to be resolved and so fluxes were extracted using 
the AIPS task JMFIT and verified using kvis. We investi¬ 
gated the flux calibration and find a typical 5% uncertainty, 
which we propagate into in our uncertainties. 

3.4 OVRO and UMRAO 

Since 2007 the Owens Valley Radio Observatory (OVRO) 
has been using its 40m telescope to undertake a 
15 GHz monitoring campaign of over 1500 radio sources 
(mainly blazar Ferm i-LAT gamma-ray candidates, see 
iRichards et al. 1 1201 ill . Five of the sources in our sample 
have been monitored as part of this OVRO monitoring cam¬ 
paign. Additionally, since January 2013 eleven extra BCGs 
with strong high radio-frequency emission selected from this 
work have been included within the dynamic queue, allow¬ 
ing regular (typically every 10 days) observations for these 
sources (see Table [I]). 

The OVRO 40 m uses off-axis dual-beam optics and 
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a cryogenic high electron mobility transistor (HEMT) low- 
noise amplifier with a 15.0 GHz centre frequency and 3 GHz 
bandwidth. The two sky beams are Dicke switched using 
the off-source beam as a reference, and the source is alter¬ 
nated between the two beams in an ON-ON fashion to re¬ 
move atmospheric and ground contamination. A noise level 
of approximately 3-4 mjy in quadrature with about 2% addi¬ 
tional uncertainty, mostly due to pointing errors, is achieved 
in a 70 second integration period. Calibration is achieved 
using a temperature-stable diode noise source to remove re¬ 
ceiver gain drifts and the flux density scale is derived from 
observations of 3C 286 assuming the iBaars et al. I (Il977i f 
value of 3.44 Jy at 15.0 GHz. The systematic uncertainty 
of about 5% in the flux density scale is not included in the 
error bars. Complete d etails of the reduction and calibration 
procedure are found in lRichards et aU\ d201lf) . We check the 
pipeline data for, and remove, individual scans whose ampli¬ 
tude is obviously discrepant (of which we find less than 1%) 
and remove a small minority of points where uncertainties 
reach more than 40%. 

Additionally, one of our sources (4C+55.16) that is in¬ 
cluded within the OVRO monitoring campaign was also 
monitored at 15 GHz by the 26m telescope of the Univer- 
sity of Michigan Radio Astronomy Observatory (UMRAO; 
lAlier et al. Ill985i. l2014h . The UMRAO monitoring ran be¬ 
tween October 1984 and June 2007, hence is not ongoing. 
However the addition of the UMRAO data to the OVRO and 
AMI data means that for this source we have well-sampled 
lightcurves at 15 GHz extending forwards from 1985 to the 
present. 


3.5 SCUBA-2 

The SCUBA-2 (iHolIand II2013T) observations were made as 
part of a poor weather programme (JCMT weather Bands 
4 and 5, tmsghz = 0.15-0.3) as part of Canadian and UK 
projects (M12AC15, M12BC18, M12BU38, M13AC16 and 
M13AU38) between February 2012 and July 2013. The ob¬ 
servations were made in small map (“CV Daisy”) mode with 
integrations of 30 minutes each. 

The observat ions were reduced using the standar d 
SMURF package lljenness et al.l[201ll : IChapin et al .11 201. 'if) . 
We used the standard flux calibration of FCFsso = 527 ± 
26 Jy beam -1 pW -1 . The resulting maps reached noise lev¬ 
els of typically 4-8 mjy beam -1 (depending on conditions) 
which was sufficient to detect over half of the sources ob¬ 
served. 

4 RESULTS 

4.1 Spectral Energy Distributions 

In H15a we performed spectral decompositions for a large 
number of BCG SEDs, finding that in many cases the spec¬ 
trum flattens above a few GHz, indicative of these sources 
containing an underlying strong core component. In the cur¬ 
rent paper we expand the spectral coverage for our subsam¬ 
ple of 35 of these sources into the mm/sub-mm regime. We 
indeed see that in many instances these active components 
extend to high frequency and typically rollover at the highest 


radio-frequencies, consistent with recently accelerated syn¬ 
chrotron populations. 

We performed fits to each of our extended SEDs, us¬ 
ing the CURVEFIT program of IDL. Generally our sources 
could be well fitted with either a power-law (equation IA2I) . 
Giga-hertz Pea ked Source (GPS) (equation IA1I also see 
lOrienti fe Dallacasa ll2014h like component or a combination 
of these. Individual SEDs, fitting notes and parameters are 
presented in Appendix m Our SCUBA-2 data were included 
on the SEDs although excluded from the spectral fits as in 
several cases it appears to be indicative of an additional 
albeit poorly constrained component becoming prominent 
in the mm/sub-mm regime (see Section [5.611 . Four example 
SEDs are shown in Figure [T) 

In some instances, as in RXJ0439+05 (panel a of Figure 
[lj the GPS-like component is distinct and inverts the spec¬ 
trum above a break frequency below which we see a steep 
spectrum power-law tail to lower frequency. These sources 
are likely to be undergoing powerful recent activity in their 
cores. At frequencies below the self-absorption turnover of 
the core component the spectrum becomes dominated by an 
underlying steep spectrum power-law component, that may 
be suggestive o f the presence of an amorphous halo of con - 
fused emission (iKempner et al. 1120 04; Hoga n et al. Il2015al l. 
In other instances, as shown in E1821+644 and 4C+55.16 
(panels b and c of Figure [T]) the integrated spectrum shows 
only a flattening in the GHz range which may stay flat to 
high frequencies as in E1821+644 or tail off as in the case of 
4C+55.16. That almost all CC-hosted BCGs contain radio- 
AGN suggests these radio-sources are long lived. The variety 
that we see in their radio spectra suggests that the emitted 
radio flux is not a constant for any given source. The timing 
of the most recent period of increased radio emission may 
then dictate the ratio of high radio-frequency core emission 
to steeper low-frequency emission from more physically ex¬ 
tended regions. Yet other sources, such as MACS0242-21 
(panel d, Figure [TJ appear well fit by a GPS model with a 
turnover frequency below 1 GHz. Such sources may be classi¬ 
fied as compact steep spectrum objects (CSS, lO’Dea Ill998lj . 
These may be dominated by a single period of older activity 
where the ejected emission has propagated outwards and ex¬ 
panded, permitting the spectral peak to move towards lower 
frequency. 

4-1.1 Notes on the General BCG Population 

Twenty-six of the thirty-five sources are found to contain 
a GPS-like peaked component. Physically, GPS sources are 
believed to be young radio sources whose spectra are peaked 
due to synchrotron self-absorption IIO’Dea Ill998h . That we 
find peaked components within our BCG spectra, in ad¬ 
dition to steep components at low-frequency, is indicative 
of recently enhanced radio emission in the central regions 
of these objects. This could be associated with a long-lived 
source having a varying accretion rate, or with the launching 
of distinct knots in the radio jets. Either scenario is consis¬ 
tent with the high duty cycle expected for these cool-core 
hosted AGN. We plot in Figure [2] both the peak flux of the 
fitted GPS-component and the combined power-law + GPS- 
component flux of each of our GPS-containing sources as a 
function of the GPS turnover frequency. These are mainly 
found to be GPS dominated at the frequency of the peak, 
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Figure 1. Example SEDs for four objects. In panel (a) we highlight the frequencies of the observations presented within this paper, 
in addition to the widely used radio-surveys NVSS/FIRST, WENSS and VLSS to help contextualise our frequency coverage. We see a 
variety of spectral shapes, sometimes well-fit by a single component as in MACS0242-21 whereas in other cases requiring both a power-law 
and a GPS-like component (see text). Note that the SCUBA-2 data-points at 353 GHz are included in the SEDs but excluded from the 
plots as they often appear to be suggesting the presence of an additional, poorly constrained component at the highest frequencies (see 
Section ESp . 


with only a minority of sources found where the combined 
flux is significantly greater than the GPS component alone. 

We find that the spectral turnover of the GPS com¬ 
ponent in nineteen of these twenty-six sources lies above 
2 GHz. Only Hydra-A has a sub-dominant core at the fre¬ 
quency of the core’s spectral peak. The presence of the core 
here therefore does not overtly affect the integrated spec¬ 
tral shape in Hydra-A at frequencies greater than a few 
GHz. In the simplest case we can use this to put a limit on 
the number of BCGs containing powerful peaked cores that 
greatly affect the spectrum above the observing frequency 
of most wide-sky surveys (e.g. NVSS/FIRST at 1.4 GHz, 
WENSS/WISH at 325 MHz, VLSS at 74 MHz, see Appendix 


Al Condon et al. 19981; 

White et al. 1997): 

Reneelink et al. 

19971: De Breuck et al. 

200^: Cohen et al. 

20Q3). Our sam- 


pie was drawn from the 726-source Parent Sample of H15a. 
If we subtract from this the 196 sources that fall below Decli¬ 
nation -30° and so were not considered by our initial GISMO 
selection criteria (see Section [5]) then we find that ^3.4% of 
BCGs contain a synchrotron self-absorbed active core com¬ 
ponent peaking above 2 GHz that is brighter than lOmJy at 


C-band (« 5 GHz). This fraction increases to ^8.5% if we 
consider only the cool-core clusters. Our incomplete spec¬ 
tral coverage of the full sample means that this number is 
only a lower limit and that the true fraction of BCGs with 
GPS-cores may be much higher. Indeed this simple analy¬ 
sis neglects systems whose spectra appear to be persistently 
flat out to high radio frequencies where multiple superposed 
flux components may be present. 

That strong spectral deviations at frequencies higher 
than a few GHz are not uncommon in BCGs has important 
implications. This includes, but is not limited to, the activity 
of the population as a whole and also the expected contam¬ 
ination rate of high-frequency peaking radio components in 
BCGs on SZ signals. We consider this further in Section [5] 


4.2 Nature of Variability 

Understanding the radio variability of AGN is important as 
it informs us on several physical processes. Short-term vari¬ 
ation allows us to place size constraints on the region from 
where any observed radio-flux propagates (i.e. a source can- 
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Figure 2. Flux of the GPS-like component as a function of the 
frequency at the peak. Note that in the classification scheme of 
O’Dea et al. (1998), the lowest peaking of these components, if 
an isolated source would likely be identified as CSS rather than 
GPS. However, we note that in the evolution scenario for this 
family of sources, CSS sources are just slightly older GPS sources 
and hence we refer to all of these peaked components as ‘GPS- 
like’ for brevity. Each source is represented twice (connected by 
dashed lines). Circles show the peak flux of the GPS component 
whereas downfacing triangles show the combined flux of both 
the GPS and any power-law component measured at the GPS 
turnover frequency. In the majority of instances we see that at 
the turnover frequency the GPS component is dominating the 
flux. It should be noted that virtually all radio-loud BCGs may 
contain a peaked active component, although in most systems 
these would be expected to peak far below the integrated flux 
and hence be undetectable in all but the brightest systems. 


• GPS Peak Flux 
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X ?s 
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not vary on less than its crossing-time) whereas longer term 
(few year) variation can inform us as to the likely fueling 
mechanisms of AGN. At radio frequencies, the variation it¬ 
self can be due to differing energy densities within the jet, 
which may be due to a change in the accretion structure and 
fuelling rate at the jet-base. Alternatively, changes in flux 
may be due to individual components within a jet becom¬ 
ing brighter as they interact with other material. However, 
fully characterising the radio variability of AGN is inher¬ 
ently difficult as the measured variability can be attributed 
to different mechanisms and will be most prevalent on dif¬ 
ferent timescales depending on the observing frequency. 

As an example, consider a synchrotron self-absorbed 
source sampled at a frequency just below the turnover which 
shows a flux increase. Without contemporary observations 
at and above the peak it would be difficult to say whether 
the additional received flux is due to an increase in total 
power-input to the jet by the AGN (and hence the normal¬ 
isation of the entire SED should increase) or whether the 
overall power-output is falling but that the lower intensity 
and expanded physical scales of the emitting regions lead 
to the turnover moving to lower frequency and so giving a 
flux increase below the peak while the overall normalisation 
remains constant. Ideally the entire SED would be sampled 
simultaneously and fitted, accounting for seif-absorption at 
a variety of epochs in order to recover the true variability 
in the underlying power output. Unfortunately, such a cam¬ 


paign is prohibitively expensive. Nevertheless, limits can be 
placed on the variability of BCGs at different frequencies 
and timescales. 


4.3 Monochromatic Variability 


4-3.1 Variability Index 

We use the de bia sed var iability equatio n (e.g . 
Akritas_&^ershadvJ Il996l : iBarvainis et al. I 120051 : 
Sadler et al. I 120061 1 to attempt to place quantitative 
limits on the monochromatic flux variability of our sources. 
This index accounts for the uncertainties in flux mea¬ 
surements and hence should be robust against artificial 
positives. The Variability Index (V.I.) is defined as 


Vrms = 


100 

<S’> 


N 


(1) 


where Si is an individual flux measurement, m the associ¬ 
ated uncertainty, < S > the mean flux and N the number of 
observations of a given source. For lightcurves with a large 
number of points this statistic gives a measure of the typical 
variability of a source about its mean flux. This analysis be¬ 
gins to falter for lightcurves with a small number of points 
where the mean is not necessarily representative of all mea¬ 
sured fluxes. However, it is useful even in the limit of N=2 
as a positive returned V.I. is indicative of significant vari¬ 
ation whereas an imaginary returned V.I. suggests any ob¬ 
served flux difference is most likely attributable to measure¬ 
ment uncertainty. A limitation of this method is that over¬ 
estimation of the error on flux measurements can mask real 
variability. However, this variability would be below any lim¬ 
its found for variation below the given errors. Given limits 
on variability would therefore be correct, albeit less restric¬ 
tive than if less conservative flux uncertainties were used. 
Although it should be noted that such a limit does not of 
course restrict the potential for the source to vary below 
this level. Similarly, under-estimation of errors can falsely 
indicate variability where there is none. 


4-3.2 Variability at 15GHz: OVRO 

For the eleven sources added to the OVRO monitoring cam¬ 
paign circa MJD 56320 (January 2013; see bottom panel, 
Table [TJ we find that no sources have a positive Variability 
Index. This suggests that for all eleven sources the difference 
in measured flux is within the l-cr uncertainties of the indi¬ 
vidual measurements such that this leads to an imaginary 
value for the square root. This is somewhat surprising, since 
for a non-varying population we might expect around half to 
have slightly positive V.I.s and half slightly negative V.I.s. 
The lack of any positive V.I.s indicates that our flux uncer¬ 
tainties are over-estimated. A positive V.I. requires that the 
numerator within the root of Equation [T] be positive, which 
at the most basic level just requires that the average frac¬ 
tional flux deviation from the mean flux be larger than the 
mean fractional error. 

In Figure [3] we plot the typical percentage deviation 
from the mean flux for each source, as well as the mean 
fractional error on the flux (displayed as a percentage). 
Additional data are provided in Table [I] We see that in 
all instances of the eleven sources monitored since January 
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Figure 5. Longer term lightcurves for the five sources that were in the original list of OVRO targets. Where AMI observations were 
available, the flux is corrected to 15 GHz using an index fitted to the AMI data split over the six channels for that observation specifically. 
We see strong variation over typically year timescales in both RXJ0439 and A2270 at 15 GHz. Note that MJD 54000 corresponds to 
September 22 nd 2006. 
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Figure 3. Dotted line shows the minimum percentage variabili¬ 
ties that would have corresponded to a positive Variability Index 
(typical deviation from mean flux greater than mean fractional 
error) for the OVRO sources, sampled at 15 GHz on typically one- 
week timescales. The green dashed line shows the actual typical 
deviation from the mean flux for each source respectively. Note 
that only A2270 shows significant variation (dashed line above 
dotted line). Sources are arranged left-to-right by their maximal 
measured flux at 15 GHz. The detectable percentage variation 
is typically much lower for brighter sources, with some devia¬ 
tion from a one-to-one correlation due to uncertainties on flux 
measurements varying by source. Sources are: A) 4C+55.16, B) 
RXJ0439*, C) A2270*, D) A2052, E) RXJ1558*, F) RXJ2341, 
G) RXJ0132, H) AS780*, I) RXJ1350*, J) A2055, K) A2415, 
L) Z8276, M) A2627, N) A2390, O) A646 and P) A1348. Note 
that five of these (indicated with *) were in the initial OVRO 
monitoring list and hence longer time-baseline data are available, 
where plotted values are those calculated when the time-frame 
was matched to the other eleven. See also Table 1T1 


Figure 4. OVRO variabilities for only the five sources in the orig¬ 
inal lists, measured for the full (from January 2008 onwards, red 
lines and offset for clarity) and restricted (January 2013 onwards, 
black lines) time-frames. Dotted lines show the minimum per¬ 
centage variabilities that would have corresponded to a positive 
Variability Index (typical deviation from mean flux greater than 
mean fractional error), whereas dashed lines show the actual typi¬ 
cal deviation from the mean flux for each source respectively. Note 
that only A2270 shows low-level significant variation (dashed line 
above dotted line) when only data during the matched time pe¬ 
riod whereas when using the full data, both A2270 and RXJ0439 
show significant variability, as would be expected by considering 
the lightcurves in Figure [5] Although not a perfect measure, this 
does show evidence that variability at 15 GHz is most prominent 
over longer time-frames, hence considering variation of BCGs on 
longer than single year timescales is important. As in Figure [3] we 
order the sources left to right by their maximal measured flux at 
15 GHz. See also Table [l] 


2013 the typical flux deviation is less than the typical error, 
hence showing why no real V.I.s were returned. For a non 
intrinsically-varying population of sources we might expect 
~68% of flux measurements to lie within one standard error 
of the mean. That our typical percentage uncertainties are 
around twice our typical flux deviation from the mean sug¬ 
gests that our errors are over-estimated by a factor of 2-3. 
We cannot steadfastly claim that our sources are intrinsi¬ 
cally non-varying and should therefore adhere to Gaussian 
statistics, hence we merely claim that they typically are re¬ 
stricted to varying by ^10% on roughly one-week timescales. 
However, we note that this limit may be a factor of 2-3 lower. 

Two of the five sources that were in the OVRO moni¬ 
toring campaign over the full period (top and middle panels, 
Table [lj) return real Variability Indices when the V.I. is cal¬ 
culated over the full timerange (A2270 and RXJ0439) but 
only one (A2270) is recovered when the V.I. is calculated 
over a timerange restricted to match that of the other eleven 
sources (see Figure a. In Figure [5] we show the 15 GHz 
lightcurves for these five sources. Clearly both A2270 and 
RXJ0439 display large flux increases over the full timescales 
considered. It is interesting to then note that the variabil¬ 
ity of RXJ0439 as shown by the increase from MJD 54500 
- 55500 would be missed if considering only the restricted 
timescale. This highlights a shortcoming of the V.I. tech¬ 


nique. If a source varies significantly (e.g. RXJ0439, Figure 
[o] MJD 54500-55500) then its average flux difference from 
the mean will be much larger than the mean flux and so 
a positive V.I. will be returned. If the source then contin¬ 
ues to be monitored but does not subsequently vary much 
(e.g. RXJ0439, MJD 55500-57000) then the average flux dif¬ 
ference from the mean will drop significantly. If data from 
monitoring over a substantially long time is combined, this 
will eventually result in the average difference from the mean 
dropping below the average error and so the measured V.I. 
is no longer truly reflective of the variability seen early on. 
The inverse of this effect can be seen by considering A2270 
in Figures [4] & [5] A2270’s lightcurve shows that it varies 
significantly. However, when only the matched timeframe 
data is used, A2270 returns a barely significant variability 
(Figure [3j. Conversely, when the longer-term monitoring is 
included then a very high level of variability is seen (Figure 
[1, which better reflects what is evident from its lightcurve.. 
These effects suggest that at 15 GHz, variability of sources 
is most likely to be seen in campaigns with high cadence 
over several year monitoring timescales. 

The measureable percentage variation that would have 
returned a real V.I. (and therefore a significant measure of 
variability above the noise) varies from source to source, 
typically being lower for brighter objects. However, by in¬ 
spection of Figure [3j we can claim that the 15 GHz vari- 
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Figure 6. Detectable variabilities that would have corresponded 
to a positive Variability Index for the GISMO sources, sampled 
at 150 GHz on approximately 6 month timescales. Sources are ar¬ 
ranged left-to-right by their maximal measured flux at 150 GHz. 
The detectable percentage variation is typically much lower for 
brighter sources, with some deviation from a one-to-one correla¬ 
tion due to uncertainties on flux measurements varying by source. 
Sources are: A) Hydra-A, B) 4C+55.16, C) A3581, D) A2270, E) 
RXJ0439, F) RXJ2341, G) A2052, H) RXJ1558, I) MACS0242, 
J) Z8276, K) A2415, L) RXJ0132, M) AS780, N) A2627, O) 
RXJ1350, P) A2055, Q) Z8193, R) A646, S) RXJ1832, T) A2390, 
U) A496, V) E1821, W) PKS0745. 


ability of these BCGs is typically restricted to ^10% on 
roughly one-week timescales (see Tabic [p. Additionally, we 
see evidence that much larger magnitude variation may be 
common i n BCGs on few-year t imescales. This mirrors the 
finding of iHovatta et al. I (l2007ll who found variability on 
few-year timescales in a sample of blazars from a greater 
than 25 year monitoring campaign with the UMRAO and 
Metsahovi telescopes. The OVRO monitoring campaign is 
ongoing hence the statistics on this longer-term variability 
will naturally improve with time. 

4-3.3 Variability at 150 GHz: GISMO 

Of the thirty-five individual sources observed with GISMO, 
twenty-three had at least one repeat observation hence al¬ 
lowing variability to be measured. Of these, and consid¬ 
ering only detections, seven sources were observed at all 
three epochs with the remaining sixteen observed twice. 
Only four sources returned real two-point Variability In¬ 
dices, namely Z8276 (V.l.=43.5 : 2 observations), A2270 (8.2 
: 3), MACS0242 (23.2 : 2) and RXJ1558 (10.2 : 3). 

However it should be noted that working out the V.l. 
for our GISMO observed sources is less informative than for 
those observed with OVRO since the V.l. works primarily 
by determining the difference from the mean flux. Clearly, 
when only two (or three) observations of a source are avail¬ 
able then such an average is less meaningful due to the large 
uncertainties on the flux. Nonetheless, a real V.l. still indi¬ 
cates which sources are really variable at the investigated 
cadence. Therefore instead of calculating the minimum de¬ 
tectable typical percentage variation about the mean (as was 


Figure 7. Ratio of the measured GISMO fluxes at 150 GHz be¬ 
tween each pair of epochs. The dashed lines indicate the error en¬ 
velope due to the 8% typical flux uncertainty. Open symbols are 
for sources that appear more than once due to being observed dur¬ 
ing all three epochs; consult also tablc[2] Sources typically appear 
to be fading, as expected for a sample chosen as the brightest at 
the selection epoch. The source that appears to show very strong 
variation is Z8276, which is seen to fade drastically at 150 GHz 
between April 2012 and April 2013. 


done for the OVRO data, see Figure [3}, for GISMO we in¬ 
stead calculate: 

i) For sources with two observations - the actual minimum 
required percentage change between the first and second 
fluxes to return a real V.l. (assuming the reported uncer¬ 
tainties) 

ii) For sources with three observations - the total minimum 
required percentage change across three epochs, from the 
first measured flux assuming a simple model where the ob¬ 
servations are symmetrically distributed around the flux of 
the central observation (i.e. Si - S 2 = S 2 - S 3 ). This is a sim¬ 
plistic approach and would require a source that does not 
vary between epochs 1 and 2 to vary by twice the expected 
mean between epochs 2 and 3 to be detected, however it 
does allow an estimate of the minimum average variation 
per six-month window that would be detectable. 

The results of this are plotted in Figure [ 6 ] We see that we 
are sensitive only to variation of above typically 30%. The 
actual percentage variations of the four sources with real 
V.L values are also shown. 

In Figure [7| we plot the fractional flux changes between 
each pair of GISMO observations for a source. As in Figure 
[ 6 ] we see that the majority of observed variability is within 
the error envelope. The large variations of some sources as 
seen in Figure[ 6 ]are again apparent. One point that becomes 
obvious in this particular representation is that the major¬ 
ity of our observed variation is of sources fading over the 1-2 
years covered. This may hint towards systematic flux offsets 
between GISMO runs. However, these are carefully checked 
foiQ and should not dominate. A natural alternative expla- 


3 Performance reports available for each run at 
http: //www. iram.es/IRAMES/mainWiki / Continuum/GISMO/M 
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Table 1. Variability Indices for the OVRO data. The mean and median timescales between observations are given in days, with N being 
the total number of times OVRO observed each source during the considered monitoring period (MJD 54470 to MJD 56750, 5 th Jan 
2008 — 3 rd Apr 2014). Note that a real V.I. is recovered for RXJ0439 when the full monitoring period is considered but not when only 
using the restricted timescale (see also Figure [5] A much better understanding of the longer-term variability of BCGs is expected to 
become apparent with continued monitoring. 


Source 

V.I. 

Corresponding % Min. detectable % Mean 

Timescale 

Median 

Timescale 

Mean 

Flux (mJy) 

N 

RXJ1558 

NaN 

Extended 

NaN 

Time Baseline (from 
5.9 

MJD 54470) 

7.5 

4.0 

195.8 

303 

RXJ0439 

8.0 

4.0 

4.0 

5.8 

4.0 

131.9 

393 

AS780 

NaN 

NaN 

10.3 

9.0 

6.0 

130.0 

205 

A2270 

22.0 

4.4 

4.2 

7.0 

4.0 

308.4 

326 

RXJ1350 

NaN 

NaN 

6.8 

7.1 

4.0 

201.2 

321 

RXJ1558 

NaN 

Matched Time Baseline (from 
NaN 5.5 

MJD 56320) 

8.3 

5.0 

261.1 

52 

RXJ0439 

NaN 

NaN 

3.8 

6.3 

5.0 

131.5 

67 

AS780 

NaN 

NaN 

10.5 

9.6 

5.9 

128.9 

45 

A2270 

1.7 

2.5 

2.5 

7.1 

5.0 

321.9 

60 

RXJ1350 

NaN 

NaN 

6.3 

8.4 

5.0 

195.5 

50 

4C+55.16 

0.48 

Added to Monitoring List MJD 56320 

2.2 2.2 6.8 

5.0 

1539.2 

63 

A1348 

NaN 

NaN 

15.9 

11.4 

7.0 

55.2 

37 

A2052 

NaN 

NaN 

6.4 

6.1 

4.2 

272.5 

71 

A2055 

NaN 

NaN 

11.8 

6.2 

4.9 

120.1 

70 

A2390 

NaN 

NaN 

11.9 

5.6 

4.8 

88.5 

76 

A2415 

NaN 

NaN 

10.6 

7.3 

5.0 

123.1 

59 

A2627 

NaN 

NaN 

11.3 

5.5 

4.2 

99.4 

77 

A646 

NaN 

NaN 

14.8 

9.1 

5.0 

58.2 

47 

RXJ0132 

NaN 

NaN 

9.3 

7.2 

5.0 

155.5 

60 

RXJ2341 

NaN 

NaN 

7.2 

6.6 

4.1 

156.9 

65 

Z8276 

NaN 

NaN 

11.6 

6.8 

4.3 

95.3 

62 


nation arises simply because we selected the sources to be 
(likely) the brightest amongst our parent sample at 150 GHz 
as a result of their lower frequency spectral shapes. These 
most active sources are therefore likely to have been selected 
at their peak brightnesses and subsequent follow-ups will 
naturally be expected to see the flux falling. 

The April 2013 observation of M0242-21 was observed 
with a high optical extinction (tau = 0.70), which may con¬ 
tribute to the GISMO flux drop of about 50%. However, we 
note that this source is a Submillimeter Array (SMA) cali¬ 
brator and shows a similar percentage flux decrease over this 
period (at 1.3nur@) suggesting that the variability is real. 

Given the high error envelopes, that we detect variabil¬ 
ity amongst our sample of repeat observations at 150 GHz 
suggests that significant variation of such sources at high ra¬ 
dio frequencies may be relatively commonplace over 6-month 
to year timescales. Such variability has important implica¬ 
tions particularly for Sunyaev-Zel'dovich observations (see 
Section [5j). 

4.4 Measuring Percentage Variation 

In addition to a V.I., we also determined a measure of the 
typical annual absolute percentage variability of the sources 
monitored by OVRO at 15 GHz. Measuring the maximal 

4 http://smal.sma.hawaii.edu/callist/callist.html 


percentage variation of a source over any given time period 
is highly susceptible to measurement uncertainty, since one 
bad unidentified outlying flux can greatly skew the result. 
Instead we effectively determine a robust measure of the 
maximal gradient of the lightcurve during the monitored pe¬ 
riod. To do this, we measure the mean of the six percentage 
differences between the 95&5th, 94&6th, 93&7th, 92&8th, 
91&9th and 90&10th percentiles (n.b. for sources only in the 
monitoring list since January 2013 with less than 100 obser¬ 
vations, we used the 5th to 95th percentile range and the 
5 unique ranges below this to a minimum separation of the 
85&15th percentiles - beyond this we automatically classi¬ 
fied any resulting measured variability as an upper-limit). 

A mean over six ranges was found, by trial and er¬ 
ror, to be a reasonable compromise between too few mea¬ 
surements being susceptible to random fluctuations but too 
many requiring the inclusion of less-separated percentile dif¬ 
ferences, over which less variation is expected even for vary¬ 
ing sources. We found that taking the 5th to 95th percentile 
range as our maximum ensured that our measurements were 
robust against outlying data-points. Additionally, averaging 
over an even number of percentile ranges allows non-varying 
sources to have variations averaging to effectively zero (note 
that in practice the probability of it actually averaging to 
precisely zero is minimal, however the probability of it ap¬ 
proaching zero is increased with an even number). 

To determine whether these measured percentage vari¬ 
ations most reliably constitute genuine variability or an 
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Figure 8. Percentage variation at 15 GHz as a function of the position of the fitted GPS peak. Note that the percentage variation has 
been normalised to the implied annual variation. This allows comparison of sources monitored over different timescales and with a variety 
of sampling timescales. Values for OVRO sources correspond to those calculated in Section [4.41 Three sources have historical observations 
from the 1970/80s at 15 GHz with the Ryle telescope (see text Section rm and we include the percentage variations over these longer 
timeframes, calculated in a similar way for comparison. Open symbols denote flux decreases over the monitored period whereas filled 
symbols denote flux increases. We see a general tendency for higher peaking sources to be more variable although this is not a strong 
trend and the ratio is highly dynamic, as highlighted by the case studies of RXJ1558-14 (shown as large crosses) and NGC1275 (shown 
as large stars - also see Section [5j. For comparison we plot as faded points the normalised percentage variations in the OVRO monitoring 
for the HFP sample of Dallacasa et al. (2000). These inhabit a similar distribution to the peaked components within our BCGs showing 
that their behaviours is not dissimilar to young radio sources elsewhere. 


upper-limit on variability, we took a measure of the two- 
point V.I. values for each of these six percentile ranges. We 
required that a minimum of four of these return real V.I. val¬ 
ues for a value to be assigned to the percentage variation, 
otherwise the measured percentage variation was taken to be 
an upper-limit. We normalise our mean measured percent¬ 
age variation values by the mean difference in years between 
the measurements used, thus recovering a measure of typical 
annual variation. Where only a limit on variation is recov¬ 
ered, we normalise by the total monitored period in years 
(roughly 1 year for the 11 sources added to OVRO list in 
January 2013). 

Three sources (4C+55.16, RXJ0439+05 and Z8193) 
have historical 15 GHz data preceding OVRO monitor¬ 
ing from either the UMRAO campaign or pointed obser¬ 
vations with the Ryle Telescope. We combine these with the 
OVRO data and find the percentage variations over longer 
timescales. 

In Figure [5] we show our calculated absolute percent¬ 
age variations at 15 GHz as a function of the peak po¬ 
sition of the fitted GPS. A mixture of sources increasing 


and decreasing in brightness is seen. A weak general trend 
for the highest peaking sources to show most variability is 
seen. Such a trend is expected, as a higher turnover fre¬ 
quency suggests self-absorption closer to the jet-base and 
hence emission from smaller scales, which can more easily 
translate to faster emitted variability. However, it must be 
noted that the GPS-peak frequency is expected to move and 
also the variability is non-constant, hence both parameters 
are expected to undergo linked evolution. We illustrate this 
dynamic relationship by including on Figure [S] evolutionary 
tracks for two sources with long-term monitoring; RXJ1558- 
14 and NGC1275. Both are seen to move extensively across 
Figure [5] but always remain within the region occupied by 
the other points. We consider these sources as case studies 
in sections IQ1 and 15.31 

For comparision with the more general galaxy popu¬ 
lation we additionally include on Figure [8] the normalised 
percentage variations at 15 GHz for the fifty sources (of 
fifty six) in the High Frequency Peaker (HFP) sample of 
iDallacasa et al. 1 (j2000i ) that are in the OVRO monitoring 
list. These show a similar weak trend to the peaked com- 
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ponents of our BCG spectra, suggesting that our peaked 
components show similar behaviour to young recently acti¬ 
vated sources. This indicates that they are associated with 
recent episodes of enhanced radio emission. 


4-4-1 A GPS link to variability? 

The behaviour of NGC1275 and RXJ1558-14 suggest that 
a direct link between GPS-peak frequency and variability is 
not present. Instead a dynamic ratio is seen between these 
parameters whereby the position (and indeed presence) of a 
self-absorption peak in relation to the observed frequency of 
a source can have a large effect on the amount of variability 
seen. 

The position of the peak appears to weakly correspond 
to how quickly a source can vary at a given frequency al¬ 
though not necessarily how quickly it does vary over any 
given epoch. Indeed, whilst the presence of a flat or peaked 
component in a BCG’s SED can be indicative of it being 
more likely to vary at high radio frequencies there does not 
appear to be a single observable proxy for the level of varia¬ 
tion. The only way to robustly remove high frequency con¬ 
taminants remains contemporaneous observation. 

It should be noted that GPS sources actually consti¬ 
tute a range of types and that many sources classified in the 
literature as GPS and HFPs are actuall y flaring flat spec¬ 
trum quasars or blazars (see Discussion bv lTorniainen et al. I 
l2007t l. The variability of these sources is naturally expected 
to be different to genuine, long-term GPS objects. Amongst 
genuine GPS sources the peaked spectral shape is main¬ 
tained over decade timescales and lon g-term variation o f 
the absolute flux level of this is seen (jAller et al. HI). 
This longer term variation (years to decades) is similar to 
what we see in the peaked components of NGC1275 and 
RXJ1558-14 and hence further supports that the variability 
herein is likely due to similar physical processes to that in 
more typical GPS objects, which may be related to opacities 
within the jet flow. 


5 DISCUSSION 

5.1 Case Study: 4C+55.16 

As mentioned in Section ECU one of our sources, 4C+55.16, 
has near-continuous 15 GHz lightcurves available for almost 
two decades, allowing us to consider its longer term activity. 
In Figure [9] we show the combined lightcurve consisting of 
UMRAO, AMI and OVRO observations. Note that our AMI 
observations have a central frequency of 16 GHz. We fit a 
single spectral index to the SED of 4C+55.16 above 10 GHz, 
recovering an index of 1.2l|3 and use this to correct our AMI 
fluxes to 15 GHz. 

A steady decline in total flux density is seen towards 
the end of the UMRAO lightcurve for this source (2007, 
MJD « 54000). Our AMI fluxes are consistent with those 
recovered with OVRO, within the errors. However, if we re¬ 
move the absolute flux calibration of the AMI observations 
then they are perhaps self-consistent with having caught the 
source as it brightens, before it is then reasonably steady 
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Figure 9. Lightcurve for 4C+55.16. UMRAO points are the 
yearly average. Note that AMI observations have a central useable 
frequency of 16.05 GHz, which was corrected to the value shown 
here using the in-band spectral index. Note that MJD 46000 cor¬ 
responds to 27th October 1984. 


throughout the OVRO monitoring. Variability of this at 
15 GHz is slow, varying by approximately 20% on decade 
timescales. This compares to significantly less variability on 
typically one-week timescales. However, from the GISMO 
observations of this source we do get faster variation at 
higher frequencies showing that although the overall SED 
may only vary slowly over several years, short-term ‘flicker¬ 
ing’ of the flux at frequencies above 100 GHz is still evident 
over much shorter periods. 


5.2 Case Study: RXJ1558-14 

RXJ1558-14 shows remarkably similar behaviour to 
NGC1275, albeit displaying around an order of magnitude 
less flux. NGC1275 is a well-studied source which constitutes 
the prototypical example of a variable BCG with a GPS- 
like component showin g stead y mm-variability over few year 
timescales (see IDutson et al. 1 120141 . for a thorough discus¬ 
sion of the variability properties of NGC1275). 

RXJ1558-14 has been monitored by OVRO since 2008 
and shows little variation at 15 GHz over this timeframe. 
However, the source has historically been used as both a 
VLA and VLBA phase calibrator and has a wealth of ob¬ 
servations available in the NRAO archive. We can therefore 
study its behaviour over long timelines. 

In the top panel of Figure [10] we show long-term 
lightcurves at both X-band (8.4 GHz) and C-band (4.8 GHz) 
for RXJ1558-14, as well as the more recent U-band 
(14.9 GHz) monitoring. Data were compiled from the lit¬ 
erature and by downloading FITS images for all available 
observations in the NRAO-archive and extracting flux mea¬ 
surements using IMFIT. We split the lightcurves into four 
time-windows, indicated on Figure 1101 a s ‘A’, ‘B’, ‘C’ and 
‘D’. In panel A fsee lWright fe Otrupcek lfl990i l we see that 
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Figure 10. Top panel: C- (~5 GHz), X- (~8) and U-Band (ps15) 
light-curves for RXJ1558-14, showing steady decline since the first 
data point, in 1972. Note that the X-band datum point in time 
window D has a border included to differentiate it from the over¬ 
lapping U-band data. Middle panel: SED where points coloured 
corresponding to the time-windows indicated in top-panel. Bot¬ 
tom panel: Illustrative fits to the GPS component of the SED for 
each of the time windows. As time goes on the GPS-peak appears 
to move both downwards and to the lower frequency. Note that 
the highest frequency points in time-window D (and to a lesser 
degree window C) appear discrepantly high. This may indicate 
a new period of activity, illustrated as a dotted orange line in 
the lower panel. Such indication suggests that a continued VLBI 
monitoring of this source may see new parsec-scale knots form 
and propagate outwards within the next few years. 


the 5-to-8 GHz spectrum was invertecjfl, although this is not 
the case at later times. 

In the middle panel we plot the SED for RXJ1558-14 
with data-points coded to correspond to these four time- 
windows. In the bottom panel we focus only on the GPS-like 
part of the spectrum, again coding the data-points although 
here we additionally show illustrative GPS-models to the 
data in each of the time windows. As we move through time 
windows A to D we see the GPS-peak appear to move both 
downwards in flux and also in turnover frequency. 

RXJ1558-14 displays double-lobed structure on parsec 
scales with an unresol ved central core, as see n at 2.224 and 
8.154 GHz in 1997 by iFev fe Charlotl (|2000l ) and more re¬ 
cently at 4.86 GHz in 2012 (see Hogan et al. 2015c). Al¬ 
though these observations are at different frequencies and 
thus cannot be used to get any secure measure on the expan¬ 
sion, it is worth noting that no new features are seen to arise 
between these observations. The lobes seen in the 4.86 GHz 
map of Hogan et al. 2015c lie approximately equidistantly 
7 milli-arcseconds either side of the central component. If 
we associate the high core fluxes observed in 1972 with the 
launching of these knots then in the intervening 40 years the 
knots travel with an apparent transverse velocity of approx¬ 
imately 1.02 c, permitting viewing angles between 45° and 
90° that are entirely consistent with the symmetry seen in 
the source. It is therefore consistent, and indeed highly likely, 
that the very high core fluxes observed in the 1970s were at¬ 
tributable to the emission of the features which are now ob¬ 
served as expanding lobes on milliarcsecond scales and the 
subsequent fall in flux is attributable to these lobes propa¬ 
gating away from the core. Such a scenario would suggest 
that variability of RXJ1558-14 (and indeed other sources) 
above a few GHz may be associated with the launching of in¬ 
dividual jet components and may precede the time at which 
these are observable given current angular resolutions. VLBI 
monitoring of self-absorbed sources that are varying above 
a few GHZ may therefore allow such flux changes to be di¬ 
rectly linked to individual parsec-scale jet com ponents, as 
has been observed elsewhere (e.g. NGC1275: ISuzuki et all 

l2012h . 

The highest frequency data in window D (and a lesser 
degree, window C) appear discrepantly high and may be 
indicative of a new period of activity that will become ap¬ 
parent at lower frequencies only as the emitting knots move 
down-jet and the self-absorption turnover peaks to lower fre¬ 
quency. VLBI monitoring of the source may therefore detect 
new knots forming and propagating outwards in the coming 
decade. 

In further analogy to NGC1275 we note that RXJ1558- 
14 has a steep spectrum power-law tail to low frequency. In 
NGC1275 this tail is associated with a 300kpc radio mini¬ 
halo. A similar structure may be present around RXJ1558- 
14. However, we note that in the TIFR GMRT Sky Sur¬ 
vey (TGSS) imaging of this source there appears to be 
reasonably symmetric structure extending almost 200 kpc 
from the core (PA « 91.8) and hence the steep spectrum 
emission may be indicative of large scale lobes from a pre- 


6 We use ‘inverted’ within this paper to refer to radio spectra with 
increasing flux to higher frequency, which is the opposite to the 
usual scenario for a typical, unobscured synchrotron spectrum. 
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vious large outburst. Interestingly, X-ray cavities are de¬ 
tecte d by Chandra in this s ource llHlavacek-Larrondo et ai] 
120131 ; iRussell (X, al 120131 1 with total physical extent of 
17.1±2.3kpc. 

Due to its relatively good long-term SED monitoring 
we consider RXJ1558-14 as a case study in Figure [5] We 
calculate the GPS-peak frequency and measured variability 
in RXJ1558-14 for each of the four epochs highlighted in 
Figure [TU] We can directly measure the 15 GHz percentage 
variability in epoch D only, during which OVRO monitor¬ 
ing data are available. For each of the other three epochs, 
we measure the C band (4.86 GHz) and X band (8.44 GHz) 
values at the epoch boundaries, using these to to extrapo¬ 
late an estimate of the 15 GHz flux. We plot the position of 
RXJ1558-14 on Figure [5] during each of these four epochs. 
Initially, in epoch A we find typically 4% variation annually 
at 15 GHz. During this epoch the peak is just below 15 GHz 
and hence the spectrum is still relatively self-absorbed at 
this frequency. In epoch B we are further in time from the 
event that caused the peaked component of RXJ1558-14’s 
spectrum. The peak has moved to lower frequency mean¬ 
ing that a steeper part of the spectrum crosses 15 GHz and 
hence contrary to the general trend we see higher variability 
during this epoch. In epochs C and D the variability con¬ 
tinues to lessen as the turnover frequency drops. It may be 
possible that a second, sub-dominant component peaked at 
much higher frequency (see Figure fTTTl bottom panel) is also 
present in the spectrum. If this component varied differently 
to a lower peaked component then it could show increasing 
flux as the other was decreasing at 15 GHz and hence could 
lessen observed monochromatic variability. 


5.3 Case Study: NGC1275 

As presented and discussed in [Dutson et al. I (l2014h . 
NGC1275 has been regularly monitored for over four decades 
and shows long term variability both in terms of its GPS- 
peak frequency and spectral normalisation. 

In order to place NGC1275 on our Figure[5]we take esti¬ 
mates of the GPS-peak frequency and 15 GHz variability in 
three epochs: 1979-83 (epoch a, Figure [8}, 1983 -2005 (epoch 
b) and 2005-2013 (epoch c) (see also Figure 5. lAbdo et alA 
1200914 . Initially, during epoch a, the peak is well above the 
15 GHz monitoring frequency and we see a sharp increase of 
around 5.5% annually. During epoch b, the self-absorption 
peak ‘rolls up’ the spectrum to lower frequency and is ac¬ 
companied by a relatively steady decline of about 3% an¬ 
nually over an extended period that may be coinciding with 
expanding features on parsec scales. Finally in epoch c there 
is another period of sharply increasing flux (almost 12% per 
year) which coincides with the peak moving to higher fre¬ 
quencies. Overall the trend is for long-term climbs followed 
by troughs, with the variability at 15 GHz dependent on 
both the position of the peak relative to this as well as the 
underlying variability in the total normalisation of the spec¬ 
trum. 


5.4 Comparisons to the general GPS population 

An important consideration to make is whether the variabil¬ 
ity and wide variety of spectral shapes that we are seeing 


in BCGs, is exclusive only to this special class of objects or 
whether they are applicable to the wider population of radio 
sources. 

There are many examples of high peaking sources in the 
literature (e.g. iRodriguez et a(.~ll2014h but large samples are 
required to fully determine how common these are amongst 
the overall radio-source population. Classifying peaked ra¬ 
dio sources is difficult. Often only non-contemporaneous ob¬ 
servations are available meaning that variability of inher¬ 
ently flat spectrum so urces can lead them to be mis-classihed 
(e.g. see discussion in lTornikoski et al. II2009I . and references 
therein). Alternatively, even when a spectrum is sampled at 
multiple frequencies and a peak observed, without follow-up 
over several years it is very difficult to determine whether 
this is a true, slowly varying GPS source or a usually flat- 
spectrum source undergoing a rapid flare. Further compli¬ 
cating the issue, it appears as if the contamination of GPS 
catalogues by variable sources and BL-LACs is dependent 
upon the GPS host galaxy; quasar-type GPS sources are 
much more commonl y mis-identifled than GPS sou rces in 
more typical galaxies dTorniainen et aZ.~ll2005l . I2007T ) . 

In the AT 2 0G su rvey of the southern sky at 20 GHz 
(|Murphv et al. 1 l201ol l , there are 3763 sources (detection 
limit of 40mJy at 20 GHz) that have simultaneous obser¬ 
vations at 5, 8 and 20 GHz. Of these, 21% are found to 
have peaked spectra, 14% show a spectral upturn over this 
range and 8% are inverted, suggesting a peak above 20 GHz. 
This fraction of inverted and peaked sources is significantly 
higher than the >3.4% of our parent sample peaking above 
2 GHz. There is however a clear selection bias towards a 
higher fraction of inverted and peaked sources in the 20 GHz 
selected sample than in our sample of BCGs detected at high 
radio frequencies but selected from a parent sample unbi¬ 
ased by radio priors. Indeed, in the 5 deg 2 sampled down to 
2.5mJy for the AT20-deep pilot survey dFranzen et al. f2014l l 
83 sources are detected and have near-simultaneous spectra 
from 1.4-20 GHz. Of these only 2.4% show a spectral up¬ 
turn, 15.7% are peaked and 6.0% are peaked above 20 GHz. 
These reduced fractions of ‘exotic’ spectral types in compar¬ 
ison to the full AT20G suggest that the area-limited nature 
of the survey coupled with the deeper detection limit means 
that more of the ‘typical’ single-spectrum sources that con¬ 
stitute the bulk of ex tragalactic radio sources are detected. 
iFranzen et al. 1 (120141 ') find that about 12% of their sources 
vary by more than 15% over 3 years at 20 GHz, which is com¬ 
parable to the levels of variability we detect for our sources 
at 15 GHz (see Figu re [8). 

iDallacasa et al. 1 (120001 ') matched 1740 sources with 
S4.9 G Hz > 300mJy in the GB6 cata l ogue (ICondon et al. I 


Il994ll to the NVSS catalogue (Condon et al. 1 119981 1. finding 
102 candidate inverted spectra. Simultaneous follow-up of 
these candidates at 1.365 - 22.46 GHz with the VLA revealed 
55 to be genuinely inverted sources whereas the remaining 47 
were flat spectrum sources whose variability meant that non- 
contemporaneous observations had caused their spectra to 
appear peaked. These 55 sources (their HFP ‘bright sample’) 
thus mean that 3.2% of their initial matched sample have 
spectral peaks >3.4 GHz. Taking the same frequency cut 
we find more than 1.9% of our parent sample peak at sim¬ 
ilarly high freque ncies. W hilst still lower than the detected 
fraction of IDallacasa et a/~l (l2000tl . if we consider that our 
sample contains BCGs detected irrespective of their radio- 
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loudness and hence contains a portion of radio-quiet objects 
then our detection samples are in reasonable agreement. We 
note that in a follow-up paper, IStanghellini et aL I (120091 1 
define an HFP ‘faint sample’ of sources with S4.9 ghz be¬ 
tween 50 and 300 mjy in the GB6 sample. Sixty-one HFPs 
with peaks at frequencies above 1.3 GHz are identified in 
this sample however they do not state the size of the parent 
sample, meaning a comparison to their detection fraction 
cannot be made. 

Overall it appears as if the peaked components of our 
BCG spectra are similar in their properties to the general 
GPS/HFP populations. However, a worthwhile point to note 
is that of the 26 sources in which we fit a GPS-like compo¬ 
nent, 20 of these (representing 76.9%) also contain a steep 
spectrum power-law to low frequencies. The usual in terpre¬ 
tatio n of peaked radio sources is that they are young dO’Dea I 
Il998l l . The presence of an accompanying steep component in 
the majority of our BCGs rules out these being truly young 
radio sources. Instead, the peaked components could be in¬ 
terpreted as due to either enhanced radio emission at recent 
times in a long-lived source, or that the radio source in BCGs 
can be restarted on short timescales. The low level variation 
that we see in our BCGs, coupled with the high duty cycle of 
detection, points to these being very long-lived radio sources 
whose radio ou tput can show quite lar ge variations in magni¬ 
tude over time. Hancock et al l : 2010 1 followed up 21 sources 
found to be inverted in the AT20G between 8 and 20 GHz 
with the Australia Telescope Compact Array (ATCA) at 40 
and 95 GHz. Of these 21 targets, 12 were found to be genuine 
peaked sources, of which 3 (25%) showed evidence of being 
restarted rather than truly young radio-sources, indicating 
them to be longer-lived sources that have undergone a recent 
episode of high activity. Although caution must be employed 
for the small numbers that we are considering, this tenta¬ 
tively supports the idea that peaked components in BCGs 
are much more commonly attributable to enhanced activity 
in long-lived sources than in the GPS population as a whole, 
which is dominated by truly ‘young’ sources. In other words, 
BCGs do not just have a single period of strong jet-emitting 
activity and then fade. They instead live for a long time, 
but remain active and effectively young at heart. 

We point out that continuous radio core activity over a 
long time is known in galaxies not selected as BCGs; for 
example, FRII radio sources, which have measured ages 
of up to 10 8 years (see e.g. iMachalsk i et al. I [2007; 
Mathews & Guo i20 1 2h - and are still being powered - 
can have cores with radio luminosities similar to those in this 
paper, with many having flat or rising spectra. However, 
if the AGN activity of BCGs has been r egulating cluster 
cores for around half the Hubble time (e.g . IVikhlinin1l2006l : 
IPratt et al. Il201(il : iMcDonald et «i.~ll2() I il l then this implies 
activity timescales over an order of magnitude longer than 
even these long-lived sources. 


5.5 Implications as an SZ contaminant 

Considering the SED shapes of our sampled sources (see 
Appendix it is clear that extrapolation of the specta 
from below 10 GHz towards the interesting range for SZ 
decrements (between roughly 15 and 200 GHz) will in many 
cases underestimate or completely overlook the contribu¬ 
tion of an active self-absorbed component. Added to this 


is the further complication that variability brings, requir¬ 
ing contemporaneous high-resolution observations to ideally 
account for contaminants (as is of course possible with in¬ 
terferometric SZ instruments such as AMI and CARMA). 

For a sample of 45 gal axy clu sters observed at 140 GHz 
with BOLOCAM, ISavers et al. I d 201 3>h used the 1-30 GHz 
spectral index in addition to limits from their 140 GHz maps 
to constrain the contamination by point sources. Although 
they concluded that typically only about 25% of the clus¬ 
ters showed a greater than 1% fractional change of the SZ 
signal, they noted that this level of contamination was much 
more prevalent in cool-core clusters (11/17, roughly 65%). 
We note that all of the BCGs in our current sample are be¬ 
lieved to lie in cool-core clusters. ISavers et al. ( 2013 1 detect 
no clear point sources at 140 GHz from their sample of 17 
cool cores, whereas from our significantly larger parent sam¬ 
ple ( 40% of our 530 cluste r parent sample are tagged as cool 
core: iHoean et al. Il2015al l we detect 32 at 150 GHz. This 
suggests ^6.0% of BCGs in all clusters exhibit bright BCG 
emission in the mm-range rising to ^15.1% if only cool-core 
clusters are considered. These values are lower-limits since 
our 150 GHz follow-up is flux limited and also incomplete in 
that some clusters are not observed at 150 GHz. The true 
level of contamination in cool-core clusters could therefore 
be even higher. We note that the prevalence of flattened or 
inverted spectral components in non cool-core clusters at fre¬ 
quenc ies below 20 GHz is much reduced fsee lHogan et al ~ I 
l2015af ). hence the level of contamination at 150 GHz is ex¬ 
pected to be significantly lower in these systems. 

Interestingly, of the ^15.1% of the CC-hosted BCGs in 
our sample, more than half (^8.5% of the total CCs) have a 
peaked component with a turnover above 2 GHz. Underesti¬ 
mating the point-source contamination will be particularly 
severe for extrapolation of the lower frequency spectrum in 
this type of object. SZ catalogues are therefore potentially 
biased against the inclusion of cool-core clusters. The true 
magnitude of line-of-sight integrated comptonization (Y sz) 
in these systems may be commonly underestimated, due to 
BCG radio emission canceling the SZ decrement. Interest¬ 
ingly this could suggest that the mass bias between X-ray 
and SZ derived cluster masses is more prevalent than previ¬ 
ously assumed. 

During their observations, the Planck Consortium 

compiled exte nsive catalogues of h i gh rad io_frequency 

sources (e.g. Planck Collaboration I l2014al . [2011b ). In 
IPlanck Collaboration I ( 2011b ! the extragalactic source 
counts agree with those o f the South Pole Telescope 
(SPT, iMocanu et al. 2QF)jb Atacama Cosmology Tele¬ 
scope (ACT, iMarriage et al. ] hoiibh and the Wilkinson 
Micr owave Anisotropy Probe (WMAP, iKomatsu et al. I 
201 lh. at the lower frequency range of Planck. However, 
Planck Collaboration ( 2011br i show that there may be a 


steepening of the typical spectral index above 70 GHz which 
could mean that contamination of the Cosmic Microwave 
Background (CMB) power spectrum by radio sources be¬ 
low the Planck detection limit may be less than currently 
thought. However, this may only apply to the bright end of 
the luminosity function. A population of GPS/HFP sources 
(in addition to a potentially ubiquitious mm/sub-mm com¬ 
ponent, see Section 15.611 in low-luminosity AGN (LLAGN) 
and BCGs could present a low level contaminant to the CMB 
power spectrum. Deep, high radio frequency radio surveys 
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are required to shed light on the luminosity function of faint 
sources in this range. 


5.6 Potential Additional sub-mm Component 

Examination of our SEDs (see Section[A]| shows that a num¬ 
ber of sources have inverted spectra between our GISMO 
observations at 150 GHz and SCUBA-2 observations at 
353 GHz. That is, the spectral slope switches from being 
falling as frequency increases at frequencies below 150 GHz 
to rising just afterwards. This is most notable in A646, 
MACS1931-26 and RXJ1504-02. The observed indices are 
too shallow (a ss 1 . 0 , in all three cases) to be so lely at¬ 
tribu table to dust-emission (expected a>~2.5, lEdge et al. I 
1199911 . However, combining the extrapolated flat spectrum 
core component with the expected dust contribution de¬ 
termined from Herschel-SPIRE observations at 250, 35 0 
and 5 00/xm (Herschel Le nsing Survey IHLSl. lEeami et al. I 
l20ld : iRawle et_aL_ ] m3) can explain the spectral upturn 
in MACS1931-26. However, there is still unaccounted for 
flux of «20 and 10 mjy at 353 GHz in A646 and R1504-02 
respectively after accounting for both these contributions, 
roughly equivalent to P 353 GHz of 1 x 10 24 W Hz -1 for 
both sources. 

An intriguing explanation to account for this ‘miss¬ 
ing flux’ is that we may be seeing an additional compo¬ 
nent in the radio/sub-mm spectrum caused by an Advec- 
tion Dominated Accretion flow (ADAF: e.g. Mahadcvan 


19971: [Naravan et al. II 19981: Ulvested &; Ho Il200ll: iDoi et al 


20051 : Naravan fe McClintock 20081 ) . These are believed 
to exist, commonly in conjunction with jets, in low ra¬ 
diative efficiency accretion sources. Multiple models exist 
that aim to explain the bolometric spectra of low lumi- 
nosit y AGN , often invoking a combination ADAF + je t 
(e.g. iFalcke I fl996l : IWu et al. 20071 : iNemmen et al. I [20141 . 
A common feature of ADAF models is the existence of a 
synchrotron self-absorbed component that typically peaks 
at mm-wavelengths. To investigate whether such a com¬ 
ponent could explain our anomalous SC UBA-2 fluxes we 
take the base model of iMahadevan I d 19971 ) for a black hole 
mass of 5 x10 9 Mq, with a viscosity parameter of 0.3 and 
a ratio of gas to total pressure of 0.5. For a ty pical B CG 
accre tion rate in the ADAF range (e.g. see iRussell et al. I 
120131 ) of 3 xlO ” 4 - 3 xlO -5 M Edd this predicts an ADAF- 
power of - 1 x lO 22 " 23 W Hz " 1 at 353 GHZ. We note that 
this an order of magnitude lower than would be inferred 
if all of the missing flux at 353 GHz in A646 and R1504- 
02 were from an ADAF component. The modeled ADAF- 
powers scale with black hole mass and accretion rate. How¬ 
ever, to align these with our missing flux would require un- 
feasibly high black hole masses or accretion rates that would 
no longer be compatible with an A DA F-like accretion struc¬ 
ture (e.g. [Done. Gierlinski fe Kubota I i2003 ). Furthermore, 
the upturned sp ectral index is expe c ted to be close to a oc 0. 4 
for an ADAF (IMahadevan 11 19971 : lAnderson et al. 1 120041) . 
Whilst our fitted core components in A646 and RXJ1504-02 
underpredict the measured 150 GHz fluxes in both instances 
(by 1.0 and 0.3mJy respectively), a simple interpolation be¬ 
tween the ‘missing fluxes’ at 150- & 353 GHz would give 
indices far too steep to be attributable to simplest case self¬ 
absorption. Alternatively, extrapolation with an index oc 0.4 


from the missing flux at 353 GHz would overpredict our mea¬ 
sured fluxes at 150 GHz. 

For the reasons given above we rule out the presence 
of a strong ADAF component to explain our missing flux. 
We note that our findings do not rule out the presence of 
an ADAF component hidden well below our observed (or 
missing) flux. The most likely explanation for the ‘missing 
flux’ then becomes a combination of low level dust emission 
and variability in the jetted core component. 


6 CONCLUSIONS 

We have expanded the radio SEDs of a sub-sample of 35 
BCGs selected from the sample of over 700 clusters con¬ 
sidered by HI5a. These sources are all hosted by cool-core 
clusters where active feedback is believed to occur, and show 
the brightest flat or inverted components in their spectra at 
high radio frequencies. These bright cores suggest that the 
central SMBH is currently accreting at an appreciable rate 
and so determining the spectra of these cores at higher fre¬ 
quency allows us to investigate possible physical origins for 
the radio emission. We considered the variability of these 
sources at both 15 & 150 GHz, enabling us to better un¬ 
derstand the amplitudes of variation in SMBHs during their 
more active phases. We find a wide variety of spectral shapes 
and in many cases see that these core components exhibit a 
spectral peak above 2 GHz, similar to the GPS/HFP popu¬ 
lation of young radio sources. Our results can be summarised 
as follows: 

• Cool-core hosted brightest cluster galaxies contain dis¬ 
tinguishable active radio core components in over 60% of 
cases (H15a). This core can become dominant at frequen¬ 
cies higher than a few GHz. We see in our current sub¬ 
sample that >15.1% of cool-core hosted BCGs contain a 
radio source greater than 3mJy at 150 GHz (equivalent to 
a radio power of « 1.2xl0 23 W Hz -1 at our median redshift 
of 0.126) and that more than half of these (>8.5% of all 
CCs) contain a distinguishable peaked component to their 
spectra with a turnover frequency above 2GHz that would 
be missed by lower frequency observations. This shows that 
the majority of CC-hosted BCGs show recent activity, of 
which a significant fraction show strong active accretion at 
any given time. 

• These core components appear to be similar to the 
GPS/HFP populations that are usually interpreted as young 
radio sources. That these peaked components are usually ac¬ 
companied by steeper spectrum power-law emission at lower 
frequencies suggests that BCGs show near-continuous activ¬ 
ity - these sources essentially enjoy a repeated youth. 

• Sources do not commonly show strong variability at 
15 GHz on short (week to month) timescales, with a typi¬ 
cal limit of < 10 % indicating relatively steady accretion over 
these timescales. We note however that there could be vari¬ 
ability that we are missing due to it being hidden within er¬ 
ror. There is an increased incidence of variation on 6 month 
timescales at 150 GHz although, accounting for measure¬ 
ment uncertainty, we find that in most cases this variability 
can only be restricted to an upper limit of <35%. That we 
do see variability on less than 6 months at 150 GHz suggests 
this emission originates within the inner 0.01 pc of the and 
hence is tracing the inner regions of the jet or accretion flow. 
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• We find that sources can show steady variation at 
15 GHz over 1-5 year timescales, although typically with 
less than 10% per annum. Variability of up to 20% pa is ob¬ 
served in a small number of sources during the most active 
periods of their lightcurves (although we note that our selec¬ 
tion of the currently brigthest systems may more favourably 
select highly varying objects). 

• This year-scale variability is found to be weakly related 
to the position of the peak in the GPS component of the 
spectrum. Whilst for any individual source the position of 
the peak at any given time does not serve as a proxy for the 
definite amount that source will vary, a higher peak does 
indicate an increased likelihood of large scale fluctuations. 

• We find that >3.4% of BCGs in our parent sample of 
530 sources at Declinations greater than 30° contain peaked 
components peaking above 2 GHz. 

• The fraction of BCGs with a peaked components peak¬ 
ing above 2 GHz increases to >8.5% if only cool-core clus¬ 
ters are considered. Overall, we find >15.1% of cool-core 
clusters contain a 150 GHz point source greater than 3mJy. 
This suggests that much more than half of all cool-core clus¬ 
ters with bright 150 GHz central point sources have spectra 
whereby even well-determined spectral indices below a few 
GHz would give very large underestimates of the flux at 
150 GHz if extrapolated and hence constitute a potential 
contaminant for SZ surveys as these sources can wipe out 
any SZ decrements if not fully accounted for. 
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Table 2: GISMO (150 GHz) fluxes of the observed clusters. 


Cluster 

Obs. Date 

RA 

DEC 

Tau 

Exposure 

(min) 

CRUSH 

Option 

Flux 

(mJy) 

FWHM 

(“) 

Peak 

(mJy) 

Map 

RMS 

Epoch 1 : April 2012 

Z235 

20/04/12 

00:43:52.30 

24:24:24.5 

0.32 

10 

-faint 

4.7 ± 1.2 

18.0 

4.4 

0.9 

RXJ0132-08 

20/04/12 

01:32:41.16 

-08:04:07.6 

0.17 

7 

-faint 

28.8 ± 4.7 

17.7 

30.9 

1.1 

MACS0242-21 

21/04/12 

02:42:36.11 

-21:32:27.3 

0.20 

2 

-faint 

35.6 ± 6.0 

18.1 

37.5 

1.8 

A496 

21/04/12 

04:33:37.82 

-13:15:46.5 

0.23 

10 

-faint 

6.1 ± 1.3 

16.4 

5.9 

0.9 

RXJ0439+05 

18/04/12 

04:39:02.32 

05:20:37.0 

0.35 

5 

-faint 

56.0 ± 9.2 

16.4 

56.5 

2.1 

PKS0745-191 

18/04/12 

07:47:31.34 

-19:17:42.7 

0.19 

25 

-faint 

4.9 ± 1.1 

21.9 

4.7 

0.8 

A646 

17/04/12 

08:22:08.86 

47:05:54.1 

0.49 

20 

-faint 

11.5 ± 2.1 

13.4 

10.9 

1.0 

4C+55.16 

17/04/12 

08:34:55.21 

55:34:20.6 

0.10 

2 

-faint 

81.7 ± 13.2 

15.6 

79.4 

1.7 

Hydra-A 

17/04/12 

09:18:05.84 

-12:05:41.1 

0.14 

4 

-faint 

266.7 ± 42. 7^. 

24.8 

179.9 

2.3 

A1348 

18/04/12 

11:41:24.26 

-12:16:37.7 

0.15 

2 

-faint 

<6.3 

- 

- 

2.1 

RXJ1347-11 

21/04/12 

13:47:31.84 

-11:45:10.9 

0.20 

20 

-deep 

<3.3 

- 

- 

1.1 

RXJ1350+09 

11/04/12 

13:50:22.01 

09:40:16.9 

0.10 

10 

-faint 

16.0 ± 2.8 

12.7 

14.9 

1.1 

A3581 

21/04/12 

14:07:30.19 

-27:01:08.6 

0.18 

2 

-faint 

71.0 ± 11.6 

16.4 

72.4 

2.1 

AS780 

18/04/12 

14:59:29.98 

-18:10:47.8 

0.10 

5 

-faint 

17.0 ± 3.2 

16.1 

16.0 

1.6 

A2052 

11/04/12 

15:16:44.50 

07:01:16.8 

0.10 

30 

-faint 

41.7 ± 6.7 

17.3 

42.5 

0.7 

A2055 

22/04/12 

15:18:46.13 

06:13:49.3 

0.26 

5 

-faint 

13.2 ± 2.5 

17.2 

12.9 

1.3 

RXJ 1558-14 

18/04/12 

15:58:22.19 

-14:10:02.7 

0.10 

5 

-faint 

39.2 ± 6.4 

14.4 

38.4 

1.5 

Z8193 

18/04/12 

17:17:19.06 

42:26:59.2 

0.13 

2 

-faint 

12.8 ± 3.1 

10.5 

13.8 

2.7 

A2270 

18/04/12 

17:27:23.71 

55:10:52.3 

0.11 

2 

-faint 

70.8 ± 4.1 

14.9 

70.0 

2.4 

Z8276 

23/04/12 

17:44:14.88 

32:59:30.8 

0.29 

5 

-faint 

35.1 ± 5.8 

17.4 

36.8 

1.3 

E1821+644 

18/04/12 

18:21:57.49 

64:20:35.6 

0.13 

10 

-faint 

5.8 ± 2.8 

13.3 

5.4 

1.4 

RXJ 1832+68 

18/04/12 

18:32:35.80 

68:48:06.3 

0.12 

10 

-faint 

10.0 ± 2.3 

11.5 

9.3 

1.7 

MACS1931-26 

21/04/12 

19:31:49.74 

-26:34:33.2 

0.12 

20 

-faint 

7.2 ± 2.1 

18.9 

7.6 

0.9 

A2390 

22/04/12 

21:53:36.74 

17:41:40.2 

0.09 

10 

-faint 

8.1 ± 1.6 

13.3 

7.7 

0.9 

A2415 

18/04/12 

22:05:38.62 

-05:35:35.1 

0.12 

5 

-faint 

31.2 ± 5.3 

19.1 

33.0 

1.9 

A2597 

18/04/12 

23:25:19.93 

-12:07:26.8 

0.11 

5 

-faint 

15.0 ± 2.8 

16.3 

15.2 

1.4 

A2627 

20/04/12 

23:36:42.29 

23:55:31.4 

0.40 

5 

-faint 

16.8 ± 3.0 

17.9 

16.8 

1.3 

RXJ2341+00 

18/04/12 

23:41:07.10 

00:18:30.9 

0.10 

5 

-faint 

41.8 ± 6.9 

16.1 

42.2 

1.5 

Epoch 2 : November 2012 









RXJ0352+19 

31/10/12 

03:52:58.97 

19:40:59.6 

0.13 

10 

-deep 

<3.6 

- 

- 

1.2 

RXJ0439+05 

31/10/12 

04:39:02.11 

05:20:42.9 

0.15 

5 

-faint 

52.0 ± 8.5 

16.6 

52.5 

1.6 

PKS0745-191 

31/10/12 

07:47:31.30 

-19:17:45.9 

0.13 

5 

-faint 

5.8 ± 1.8 

17.8 

5.3 

1.5 

A646 

31/10/12 

08:22:10.05 

47:05:54.9 

0.10 

6 

-faint 

11.4 ± 2.1 

18.4 

12.0 

1.0 

4C+55.16 

31/10/12 

08:34:55.26 

55:34:22.9 

0.13 

5 

-faint 

77.8 ± 12.5 

18.0 

80.1 

1.4 

Hydra-A 

31/10/12 

09:18:05.83 

-12:05:46.5 

0.11 

5 

-faint 

273.8 ± 43.8 7 

24.5 

167.8 

1.8 

A1084 

31/10/12 

10:44:33.11 

-07:04:11.1 

0.11 

10 

-faint 

4.5 ± 1.2 

19.0 

4.1 

1.0 

A1348 

31/10/12 

11:41:24.37 

-12:16:40.0 

0.11 

5 

-faint 

6.9 ± 1.7 

10.7 

6.5 

1.3 

A1774 

31/10/12 

13:41:05.19 

39:59:47.1 

0.15 

10 

-faint 

6.9 ± 1.5 

19.7 

6.9 

1.0 

RXJ1347-11 

31/10/12 

13:47:31.84 

-11:45:10.9 

0.12 

10 

-deep 

<3.0 

- 

- 

1.0 

A1795 

31/10/12 

13:48:52.25 

26:35:32.2 

0.19 

10 

-faint 

5.1 ± 1.4 

17.7 

5.4 

1.1 

RXJ 1350+09 

31/10/12 

13:50:22.04 

09:40:07.9 

0.19 

10 

-faint 

16.0 ± 2.7 

17.0 

15.7 

1.0 

A3581 

31/10/12 

14:07:30.00 

-27:01:07.8 

0.12 

5 

-faint 

62.1 ± 10.1 

19.0 

63.2 

2.0 

A1885 

31/10/12 

14:13:43.86 

43:39:40.0 

0.10 

10 

-deep 

4.0 ± 1.0 

12.4 

4.2 

0.8 

AS780 

31/10/12 

14:59:29.14 

-18:10:51.3 

0.11 

10 

-faint 

15.4 ± 2.7 

17.9 

15.3 

1.2 

A2052 

31/10/12 

15:16:44.64 

07:01:16.3 

0.11 

5 

-faint 

38.0 ± 6.3 

16.3 

38.4 

1.7 

A2055 

02/11/12 

15:18:46.19 

06:13:56.5 

0.50 

10 

-faint 

10.4 ± 2.5 

18.9 

10.2 

1.8 

RXJ 1558-14 

31/10/12 

15:58:22.22 

-14:10:03.1 

0.11 

5 

-faint 

27.8 ± 4.8 

15.1 

27.6 

1.9 

NGC6338 

31/10/12 

17:15:23.25 

57:24:41.5 

0.12 

10 

-faint 

7.3 ± 1.5 

11.0 

OO 

1.0 

Z8193 

31/10/12 

17:17:19.11 

42:26:59.2 

0.11 

10 

-faint 

15.6 ± 2.8 

15.9 

15.0 

1.2 

A2270 

31/10/12 

17:27:23.52 

55:10:52.5 

0.11 

5 

-faint 

56.2 ± 9.1 

18.0 

58.3 

1.3 

Z8276 

02/11/12 

17:44:14.82 

32:59:27.4 

0.53 

10 

-faint 

12.5 ± 2.7 

16.1 

13.2 

1.8 

E1821+644 

31/10/12 

18:21:57.00 

64:20:32.9 

0.11 

10 

-faint 

7.9 ± 1.7 

11.3 

7.8 

1.2 


7 Resolved source, integrated flux 
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Table 2 - Continued 


Cluster 

Obs. Date 

RA 

DEC 

Tau 

Exposure 

(min) 

CRUSH 

Option 

Flux 

(mJy) 

FWHM 

(“) 

Peak 

(mJy) 

Map 

RMS 

RXJ1832+68 

31/10/12 

18:32:35.08 

68:48:02.2 

0.12 

7 

-faint 

8.9 ± 1.7 

17.2 

9.0 

1.0 

Epoch 3 : April 2013 










RXJ0132-08 

11/04/13 

01:32:41.05 

-08:04:07.8 

0.83 

5 

-faint 

27.6 ± 13.7 

17.4 

24.0 

13.0 

M0242-21 

11/04/13 

02:42:35.99 

-21:32:24.9 

0.70 

10 

default 

17.9 ± 6.6 

15.4 

22.9 

6.0 

A496 

10/04/13 

04:33:37.75 

-13:15:40.7 

0.57 

20 

-faint 

4.4 ± 1.1 

13.1 

4.3 

0.8 

RXJ0439+05 

10/04/13 

04:39:02.31 

05:20:42.7 

0.69 

5 

-faint 

47.6 ± 7.8 

17.4 

48.6 

1.6 

A646 

10/04/13 

08:22:09.50 

47:05:53.4 

0.59 

10 

-faint 

12.2 ± 2.3 

10.7 

11.4 

1.3 

4C+55.16 

10/04/13 

08:34:54.79 

55:34:19.8 

0.63 

5 

-faint 

61.7 ± 10.0 

16.9 

62.7 

1.8 

Hydra-A 

10/04/13 

09:18:05.79 

-12:05:39.6 

0.61 

5 

-extended 

201.9 ± 32.9 7 

25.0 

153.5 

6.2 

RXJ 1504-02 

11/04/13 

15:04:07.43 

-02:48:15.8 

0.59 

10 

-faint 

5.5 ± 1.4 

11.2 

5.1 

1.1 

A2052 

11/04/13 

15:16:44.48 

07:01:18.9 

0.58 

5 

-faint 

35.0 ± 5.8 

17.0 

36.6 

1.6 

RXJ 1558-14 

11/04/13 

15:58:22.01 

-14:09:58.2 

0.63 

5 

-faint 

25.1 ± 4.3 

14.9 

23.8 

1.6 

A2270 

11/04/13 

17:27:23.58 

55:10:52.8 

0.64 

5 

-faint 

50.4 ± 8.2 

14.5 

51.2 

1.7 

Z8276 

11/04/13 

17:44:14.47 

32:59:29.4 

0.68 

5 

default 

<44.7 

- 

- 

14.9 

A2390 

11/04/13 

21:53:36.96 

17:41:40.4 

0.65 

10 

-faint 

8.6 ± 2.0 

14.1 

8.6 

1.4 

A2415 

11/04/13 

22:05:38.64 

-05:35:31.7 

0.65 

5 

-faint 

26.7 ± 4.7 

16.3 

25.8 

1.9 

A2627 

11/04/13 

23:36:41.84 

23:55:28.2 

0.66 

10 

-faint 

16.9 ± 3.0 

14.9 

15.9 

1.3 

RXJ2341+00 

15/04/13 

23:41:07.01 

00:18:33.0 

0.39 

5 

-faint 

33.0 ± 2.9 

15.2 

32.5 

1.8 
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Table 3: AMI (16 GHz) fluxes of the observed BCGs. 


Cluster 

Obs. Date 

Flux 

RMS 

Obs. Date 

Flux 

RMS 

Obs. Date 

Flux 

RMS 

Z235 

11/04/12 

17.7 ± 0.9 

0.2 

04/06/12 

16.7 ± 0.9 

0.3 

- 

- 

- 

RXJ0132-08 

17/04/12 

139.5 ± 8.3 

4.5 

02/06/12 

126.1 ± 6.9 

2.7 

- 

- 

- 

RXJ0439+05 

12/04/12 

324.8 ± 16.3 

1.6 

30/05/12 

307.5 ± 15.5 

2.1 

29/09/12 

299.9 ± 15.1 

2.0 

A646 

10/04/12 

45.6 ± 2.3 

0.5 

30/05/12 

45.7 ± 2.7 

1.4 

27/09/12 

47.3 ± 2.5 

0.7 

4C+55.16 

14/04/12 

1317.7 ± 66.0 

3.5 

30/05/12 

1372.3 ± 70.6 

16.8 

- 

- 

- 

RXJ1350+09 

24/04/12 

133.7 ± 6.7 

0.9 

14/06/12 

107.4 ± 5.6 

1.5 

- 

- 

- 

A2052 

24/04/12 

231.8 ± 11.8 

2.2 

14/06/12 

220.8 ± 11.2 

1.9 

- 

- 

- 

A2055 

24/04/12 

66.7 ± 3.5 

1.2 

14/06/12 

68.6 ± 3.6 

1.2 

- 

- 

- 

Z8193 

24/04/12 

79.8 ±4.1 

1.1 

26/06/12 

79.3 ± 4.1 

1.2 

- 

- 

- 

A2270 

24/04/12 

185.0 ± 9.4 

1.6 

26/06/12 

253.8 ± 12.8 

1.4 

26/09/12 

226.9 ± 11.6 

2.6 

Z8276 

14/04/12 

88.6 ± 4.5 

0.7 

31/05/12 

95.9 ± 4.9 

0.8 

26/09/12 

87.6 ± 4.6 

1.3 

E1821+644 

25/04/12 

9.0 ± 0.5 

0.2 

30/06/12 

10.8 ± 0.8 

0.2 

- 

- 

- 

RXJ1832+68 

19/05/11 

45.6 ± 2.3 

0.4 

25/04/12 

42.2 ± 2.1 

0.1 

05/07/12 

42.3 ± 2.1 

0.3 

A2390 

14/04/12 

68.5 ± 12.4 

0.8 

31/05/12 

66.4 ± 3.4 

0.9 

- 

- 

- 

A2415 

11/04/12 

77.4 ± 4.2 

1.7 

02/06/12 

77.2 ± 4.5 

2.3 

- 

- 

- 

A2627 

14/04/12 

83.1 ± 4.2 

0.7 

31/05/12 

73.3 ± 4.0 

1.5 

- 

- 

- 

RXJ2341+00 

11/04/12 

147.0 ± 7.4 

1.2 

31/05/12 

140.9 ± 7.2 

1.5 

28/09/12 

209.2 ± 11.8 

5.4 
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Table 4: CARMA (90 GHz) fluxes of the observed BCGs. 


Cluster 

Obs. Date 

RA 

DEC 

Flux 

RMS 

Z235 

11/06/12 

00:43:52.21 

24:24:21.5 

9.9 ± 1.8 

4.5 

RXJ0132-08 

21/05/12 

01:32:41.11 

-08:04:05.7 

98.8 ± 11.9 

6.7 

MACS0429-02 

15/06/12 

04:29:36.00 

-02:53:07.0 

<14.1 

4.7 

A496 

15/06/12 

04:33:37.89 

-13:15:42.3 

11.4 ± 3.6 

3.4 

RXJ1350+09 

21/05/12 

13:50:22.12 

09:40:10.7 

22.8 ± 2.8 

1.7 

A3581 

21/05/12 

14:07:29.83 

-27:01:04.2 

82.5 ± 9.1 

2.8 

AS780 

21/05/12 

14:59:28.78 

-18:10:44.9 

29.5 ± 4.4 

2.3 

A2052 

21/05/12 

15:16:44.58 

07:01:18.1 

126.7 ± 14.3 

9.1 

A2055 

21/05/12 

15:18:46.45 

06:13:57.9 

53.6 ± 9.0 

21.6 

RXJ1504-02 

21/05/12 

15:04:07.50 

-02:48:16.9 

8.3 ± 1.6 

1.4 

RXJ1558-14 

21/05/12 

15:58:21.91 

-14:09:58.9 

44.8 ± 4.9 

6.3 

Z8193 

22/05/12 

17:17:19.21 

42:26:57.8 

22.4 ± 2.9 

1.8 

A2270 

22/05/12 

17:27:23.49 

55:10:53.9 

164.7 ± 17.2 

14.4 

Z8276 

22/05/12 

17:44:14.47 

32:59:27.4 

35.4 ± 6.0 

4.9 

E1821+644 

24/05/12 

18:21:56.24 

64:20:58.0 

7.2 ± 1.8 

1.7 

RXJ1832+68 

22/05/12 

18:32:35.52 

68:48:07.2 

19.6 ± 4.0 

3.5 

MACS1931-26 

11/06/12 

19:31:49.64 

-26:34:32.4 

9.6 ± 2.9 

2.7 

A2390 

04/06/12 

21:53:36.81 

17:41:44.8 

22.3 ± 3.2 

2.3 

A2415 

04/06/12 

22:05:38.53 

-05:35:33.7 

25.7 ± 3.4 

2.3 

A2597 

21/05/12 

23:25:19.82 

-12:07:28.6 

19.7 ± 5.8 

2.0 

A2627 

04/06/12 

23:36:42.82 

23:55:23.8 

19.8 ± 4.5 

4.0 

RXJ2341+00 

21/05/12 

23:41:06.94 

00:18:33.3 

72.8 ± 8.0 

3.2 
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Table 5: SCUBA-2 (353 GHz) fluxes of the observed BCGs. 
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Cluster 

Obs. Date 

Flux (mjy) 

RMS 

RXJ0132-08 

23/10/12 

- 

13.1 

MACS0242-21 

01/10/12 

34.0 

7.6 

RXJ0439+05 

01/10/12 

18.7 

5.2 

A646 

27/02/12 

25.9 

6.6 

4C+55.16 

27/02/12 

29.0 

6.5 

Hydra-A 

08/10/12 

76.3 

8.8 

A1348 

27/12/12 

- 

7.9 

RXJ1347-11 

29/12/12 

- 

5.5 

A1795 

27/12/12 

- 

7.4 

RXJ1350+09 

23/02/12 

- 

5.7 

A3581 

28/01/13 

59.2 

9.1 

AS780 

15/04/12 

- 

5.0 

RXJ1504-02 

12/01/13 

13.7 

4.0 

A2052 

28/01/13 

- 

7.8 

RXJ1558-14 

15/04/12 

40.3 

4.7 

Z8193 

14/05/12 

- 

8.5 

A2270 

30/09/12 

34.9 

4.7 

Z8276 

30/09/12 

27.7 

5.0 

RXJ1832+68 

05/09/12 

- 

5.8 

MACS1931-26 

16/09/12 

16.8 

4.8 

A2390 

24/04/13 

- 

4.9 

A2415 

01/10/12 

18.2 

6.7 

A2597 

24/10/12 

- 

7.2 

A2627 

24/10/12 

20.1 

7.3 

RXJ2341+00 

31/08/12 

33.2 

3.9 
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APPENDIX A: NOTES ON INDIVIDUAL 
SOURCES 

Individual source SEDs were populated using the main 
radio catalogue s (e.g. Australia Tele scope 20 GHz Sur¬ 
vey (AT20G), iMurphv et al. I l201(j : NVSS and Faint 
Images of the Radio Sky a t Twenty-cm (FI RST) at 


1.4 GHz, Condon et al. 199£ 

; White et al. 1993; SUMSS 

at 843 MHz. [Mauch et al. 

2003j; Texas Survey of Ra- 


dio Sources (TEXAS) at 365 MHz, iDouglas et al. 11 19961 : 
Westerbork Northern Sky Survey (WENSS) and Wester- 
bork In the Southern Hemisphere (WISH) a t 325 MHz, 
iRengelink et al ] Il997l : iDe Breuck et aC\ 120021: VLA Low- 
Frequ ency Sky Survey (VLSS) at 74 MHz. Icohen et aL I 
[ 2 OO 3 ), the NASA/IPAC Extragalactic Database (NED) and 
the High Energy Astrophysics Science Archive Research 
Center (HEASARC) database, the National Radio Astron¬ 
omy Archive (NRAO), literature searches, the data pre¬ 
sented herein and supplemented by data from H15a. 

Due to the lack of self-similarity between these spectra 
each source was fit on an individual basis. Where possible, 
act ive components were isolate d and fit using the GPS model 
of dOrienti fc Dallacasa II 20 T 4 I . :) 

Log(S) = Ao + Log(y)(A\ + A 2 Logr(^)) (Al) 

Steep spectrum components and active flat spectrum com¬ 
ponents where no peak was apparent were both fit using a 
simple power-law of the form: 

S = A v~ a (A2) 

Where both an active and inactive component were found 
to be present, the core contribution is removed from the 
total emission prior to fitting the more extended emission. 
In a couple of cases clear spectral curvature attributable 
to synchrotron ageing is present in the steeper spectrum 
components and is modelled as described in the notes below. 
Additionally, flattening of a steep spectrum component 
to low frequencies (below 100MHz) is seen in a minority 
of sources. Where morphologically it is clear that this is 
extended emission, and hence the flattening is most likely 
due to free-free absorption then these data are excised from 
the fits as described below. 

4C+55.16 

2MFGC 06756, z=0.2411. This is a powerful radio source 
that is clearl y extended on mil liarcsecond scales as observed 
using VLBI (iHelmboldt II2007I . see also Hogan et al. 2015c, 
in prep.). The SED has a flat index in the GHz range 
that appears to be caused by the superposition of a broad 
peaked GPS-like component with a steep spectrum power- 
law component to lower frequency (below about 100 MHz). 
The SED is fit using a GPS-like model with parameters 
A 0 =3.90±0.03, Ai=-0.04±0.02 and A 2 =-0.43±0.01 and a 
steep power-law to lower frequency with a=l.l±0.1 and 
A=500.0±10.0. 

A1084 

2MASX J10443287-0704074, z=0.1326. The SED is sparsely 
samp led, although the non-detection in VLSS J Cohen et al. I 
120071 1 coupled to its flat returned index shows the source 
to be core-dominated. The SED is fit with a power-law of 
a=0.36±0.04 and A=36.7±2.9. 

A1348 


LCRS B113851.7-115959, z= 0.1195. The source is best fit 
using a GPS-like model, with parameters Ao = 1.96±0.20, 
Ai=0.62±0.41 and A 2 =-0.534=0.20 and then a steeper spec¬ 
trum component that becomes apparent below 1 GHz that 
is fit using a power-law with a=l.054=0.07 and A=80.24=1.5. 
A1774 

2MASX J13410515+3959456, z=0.1715. The SED of this 
source is consistent with a single power-law fit with param¬ 
eters of a=0.614=0.03 and A=104.94=5.0 , that persi sts to 
high frequency. Note that the FIRST (I White et al. 1 111)971 ) 
datum-point was removed from this fit as it appears to have 
resolved out some emission. The source appears one-sided 
in FIRST but the core is not isolated and the lack of 
higher-resolution data therefore means that a separate core 
component cannot be fit. The intermediate spectral index 
combined with the resolved morphology suggests that an 
active core is likely to be present in this system, which may 
be distinguishable with higher resolution observations. 
A1795 

CGCG 162-010, z=0.0632. Extended source with sev¬ 
eral observations of high enough resolution to cre- 
ate a separate SED for the i sol ated core on ly (see 
Laurent-Muehleisen et al. ] 1 19971 : iLin fe Moh r 120091 : 

Liuzzo et al. I 120091 b This core spectrum is flat to high 
frequency, and fit using a power-law with a=0.654=0.06 and 
A=120.2±14.2. The SED of the more extended emission 
appears to exhibit spectral curvature to high frequency, 
hence to account for this we allow an exponential rollover 
to the steep power-law, of the form: 


S = Au~ a (l- Be~^) (A3) 

llHogan et al. I l2015al h The best fit returns parameters 
A=1262.8, a=0.86, B=1.07 and ^0=7.51 for the unknown 
parameters in this model. 


A1885 

2MASX J14134379+4339450, z=0.090. The SED is core 
dominated with a hint of vari a bility seen from hi st orical 
C-band observations (lLaurent-Muehleisen et al. I Il997l : 
lGregorv_j£_Condoir _ 199 ill . Additionally, the WENSS 


( Rengehnf^etaL 


1997l l detection alongside the VLSS 


( Cohen et al. 2007) non-detection shows that the spectrum 


must peak in the GHz range. The SED can be fit with a 
GPS model, returning parameters of 1.66±0.03, 0.24±0.04 
and -0.33±0.03. 

A2052 


UGC 09799/3C 317, z=0.0355. A separate core only 
SED can be created for the resolved isolated core. This 
can be fit with a GPS model (returning parameters; 
2.54=1=0.02, 0.104=0.02 and -0.24±0.01). Accounting for this 
core contribution, fit a power-law (with a=1.20±0.08 and 
A=6404.54=1142.8) to the lower frequency emission. Note 
that there are a number of very low frequency (less than 
50MHz) data-points that are flatter than this steep spec¬ 
trum would imply. These are possibly affected by free-free 
absorption and hence were excised from the power-law fit. 
A2055 

2MASX J15184574+0613554, z=0.1019. High resolution 
data reveals an underlying core component with a GPS-like 
spectrum. Note that the CARMA datum is tagged as 
potentially unreliable and indeed appears inconsistent and 
has therefore been excluded from the fit. A GPS fit to the 
core-only SED returns parameters of 1.28±0.04, 1.134=0.02 
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and -0.57+0.01. The integrated SED is increasingly domi¬ 
nated by a power-law component at frequencies less than 
1 GHz. Fit this separately, accounting for the expected core 
contribution, to recover a power-law fit with a=0.70±0.04 
and A=519.2±122.7. 

A2270 

2MASX J17272346+5510538, z=0.2473. Source appears 
highly vari able. There is a l imit only in the VLSS-redux 
catalogue (lLane et al. [Hi) but an overlay of the map 
shows there to be a 3<r map detection at the optical position 
of the B CG. Along with other low frequ ency detections 
from 7C (lHales et al. Il200i 'l and WENSS dRengelink et al~ I 
H), this hints at there being a low-powered low frequency 
power-law component within this source. The data are 
widely spaced in observation date, making separate fits 
by epoch for the GPS component unfeasible. Fitting a 
GPS to fit to all the available data at higher frequencies 
than L-band, recovers an average GPS with parameters 
A 0 =2.11±0.06, Ai=0.56±0.06 and A 2 =-0.31+0.04. Ac¬ 
count for the contribution of this to the lower frequency 
data and find a resultant power-law component with 
a=0.94±0.07 and A=28.0±0.7 
A2390 

ABELL 2390:[YEA96] 101084, z=0.2328. The spectrum 
appears to have a split power-law like shape below 10 GHz 
although the higher frequency observations show that the 
active component curves. Hence the SED is best fit using a 
combination of a GPS component (that returns parameters; 
2.28±0.03, 0.26±0.04 and -0.42±0.05) and then a power- 
law compon ent with a=1.34±0.25 and A=121.5±60.6. 
iSmail et al. I (j2002l f claim a 6.7mJy SCUBA detection at 
330 GHz for this source. This suggests that there may be 
additional intermittent activity at these higher frequencies 
or alternatively there may be a dust component coming 
in (although non-detection (<60mJy) at 660 GHz opposes 
this). Note that we do not include this claimed detection 
on our SED, including only our more contemporaneous 
upper-limit at 330 GHz. 

A2415 

2MASX J22053865-0535330, z=0.0573. The SED has 
a flat index from around 1 GHz up to high frequency 
although a subtle drop-off means that this is best fit with a 
broad-peaked GPS model (returned parameters: 2.09±0.08, 
0.26±0.02 and -0.29±0.01) in addition to a steep power-law 
component at lower frequencies with a=1.31±0.02 and 
A=193.2+9.4. 

A2597 

PKS 2322-12, z=0.0830. Source appears very flat to higher 
frequencies. Archival VLA observations at C, X and 
U-bands from the NR.AO Image Archive allow isolated 
core measurements which in conjunction with higher 
frequency data suggest a broad GPS-like core that is well 
fit with parameters Ao=1.92±0.05, Ai=0.27±0.09 and 
A 2 =-0.30±0.04, although there may be further flattening 
to yet higher frequencies as suggested by SCUBA observa¬ 
tions (IZemcov et al. I 120071 : iKnudsen et al. I 120081 '). Lower 
frequency data are fit with a power-law of a=l.18+0.06 
and A=2336.4±13.4, where the data below 100 MHz have 
been excluded as here the spectrum flattens, possibly due 
to free-free affected. Note that a simple extrapolation of 
this power-law in conjunction with the GPS core gives a 
wholly inaccurate integrated spectrum and hence the older, 


extended component must roll-off at some frequency. This 
has been simulated here by subtracting the GPS from the 
simple power-law to give the presented curved spectrum. 

A2627 

This cluster has two bright galaxies in close proximity, 
both of which are radio-sources. The more dominant source 
- A2627(a) / B2 2334+23 ‘the BCG’ - is at position 
23:36:42.12 23:55:29.0 and z=0.126. The other galaxy - 
A2627(b) / 2MASX J23364245+2354442 is at position 
23:36:42.47 23:54:44.6 and has a redshift given in NED as 
z=0.122. VLA C-band data (H15a) separates the sources 
and overla ying other maps ont o this we see that at NVSS 
(1.4 GHz. lCondon et al. 11 1 9981 ) and lower frequencies the 
sources are confused. They are similarly confused in our 
AMI maps. The BCG can be isolated in the CAR.MA map 
and a VLA-A array image from the NR.AO Image Archive. 
These show it to be flat spectum. Along with the GISMO 
data, fit a flat spectrum power-law to this (a=0.37+0.06 
and A=109.0±22.9). Morphologically, A2627(a) is core 
dominated although it does show some wispy emission that 
is likely to contribute at lower frequencies. A2627(b) is 
fainter (by roughly a factor of 4 at C-band) but appears 
to be steep spectrum and is likely to become increasingly 
dominant at lower frequency. Fit the integrated emission 
(accounting for the flat component) with a power-law of 
a=0.65±0.06 and A=499.8±6.4 noting that this is highly 
confused. We currently do not have sufficient coverage to 
fully deblend the sources. 

A3581 

IC 4374, z=0.0218. Flat spectrum core dominated system 
with hints of variation in repeat observations around 
C- and X-bands. Unresolved at VLA-A in a map from 
the NR.AO Image Archive but the source is resolved on 
milli-arcsecond scales with the VLBA (Hogan et al. 2015c, 
in prep.). This small-scale resolved emission is linked to 
ongoing activity and the integrated flux recovered with the 
VLBA is consistent with lower resolution fluxes showing 
that source is core-only. Fit a power-law to the integrated 
SED (returning a fit of a=0.49±0.04 and A=758.2±6.9) 
and then highlight the VLBA peaks on the SED. 

A496 

MCG -02-12-039, z=0.0328. This source has a GPS-like 
core with a steeper power-law component to low frequency. 
The core is fit with the GPS-like model with parameters 
A 0 =1.45±0.02, Ai=0.60±0.02 and A 2 =-0.43±0.03 and 
then accounting for this component, the remaining ex¬ 
tended emission is fit with a power-law of a=1.71±0.03 and 
A=99.8+6.4. 

A646 

2MASX J08220955+4705529, z=0.1268. The SED of this 
source can be fit using a GPS model for the core compo¬ 
nent, which returns parameters 1.89+0.02, 0.14+0.01 and 
-0.25+0.08 and a power-law component of a=1.54±0.09 
and A=13.5±0.3 to account for the lower frequency, 
extended emission. Note that the combination of these 
two components at 1 .4 GHz over-estimates the NVSS flux 
dCondon et al. |[l998li which may be attributable either to 
variation of the core component or perhaps indicating that 
this source should in reality have a more sharply peaked 
profile. 

AS 780 

2MASX J14592875-1810453, z=0.2344. This object has 
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previously be en id entified as a flat spectrum source (e.g. 
iHealev et al. Il2007l ). The GHz range flat spectrum falls off 
above about 20 GH z. The source is undetected in VLSS 
dCohen et al. I 120071) . which suggests that the spectrum 
must turnover around 1 GHz. The SED is fit with a GPS 
model, that returns parameters 2.16±0.02, 0.34±0.02 and 
-0.35T0.01. Note that the source is m ildly confused in the 
WISH map (iDe Breuck et al 200l). The flux used here 
is deblended BCG emission however there could still be 
wrongfully attributed flux contributing to this 352 MHz 
measurement. Equally, this 352 MHz point lying above the 
curve could be indicative of a steeper spectrum component 
(which would push the GPS peak to higher frequency). 
Future low frequency observations will determine this. 
E1821+644 

HB89 1821i643, z=0.2970. This source is one of the few 
radiatively effic ient AGN k no wn to reside within a cool core 
cluster (see e.g. iRussell et 0^71120101 . and references therein). 
Nevertheless, the radio-properties appear to be consistent 
with its more inefficent brethren. Multiple high resolution 
VLA observations from the NRAO Image Archive allow 
the core to be isolated at multiple frequencies which are 
consistent with the spectral flattening seen above about 
10 GHz. The source shows significant variation. Fit the 
core-only measurements with a GPS-like model, returning 
parameters of 0.83T0.10, 0.40±0.06 and -0.19±0.03. This 
suggests a broad, flat profiled core. Independently fit a 
steep power-law to the low frequency data, returning a fit 
with a=1.06±0.13 and A=94.6±3.6. 

Hydra-A 

This source is associated with the bright radio source 
3C 218, z=0.0549. This source is very we ll studied (e.g. 
iMcNamara et qZ~l l2000l : lHamer et al. 1 120141 . and references 
therein). VLB A observations at multiple frequencies (e.g. 
iTavlor I|l996l : lArava et al. II 2 OI 0 I ) allow the SED of the core 
only to be independently fit with a GPS model (returning 
parameters: 2.01±0.05, 0.95±0.06 and -0.59±0.06). The 
integrated SED is dominated by flux from the extended 
emission and remains persistent to frequencies higher than 
100 GHz. Accounting for the minimal core contribution, 
fit the SED using a power-law with a=1.00±0.02 and 
A=47214.3±43.0. 

MACS0242-21 

PKS 0240-217, z=0.314. This source appears to be a CSS 
or low-peaking GPS, that can be fit using a GPS model 
with parameters of Ao=3.14±0.02, Ai=-0.14±0.02 and 
A2=-0.3 0±0.02. VLBA observ ations of the source at S- and 
X-band (jBeaslev et al. I l2002i ) show that it is ‘wispy’ on 
parsec scales, with the integrated flux consistent with the 
unresolved emission at VLA resolutions. That the source 
has a low spectral peak, alongside it being very structured 
on parsec scales, is consistent with this source being in a 
late-stage of core activity. 

MACS 1931-26 

PMN J1931-2635, z=0.3520. There are two radio-sources 
in close proximity within the cluster core; the BCG at 
19:31:49.67 -26:34:33.4 and a Wide-Angled Tail (WAT) 
source, PMN 1931-2635, at 19:31:50.0 -26:35:16.4. These 
sources are confused in many radio maps of the region. 
The BCG core is isolated in ATCA-6km observations at 
5.5 and 9.0 GHz along with CARMA and GISMO at 
higher frequencies. The SED of only these points is well 


fit with a GPS model, returning parameters 0.95±0.04, 
0.42±0.03 and -0.21±0.01. The sources are confused in 
lower resolution observations. The sources cannot therefore 
be fully distinguished, however by fitting to the flux peaks 
coincident with the optical positions it appears as if PMN 
1931-2635 has a steeper spectrum (as would be expected 
for a lobe-dominated WAT). Further confounding the 
picture, hig h resolu tion VLA-A observations at L-band by 
lEhlert et al. I (1201 ll ) show that the BCG is surrounded by 
an amorphous halo (see also iMittal et al. 1 120091 ) that is 
likely to be very steep spectrum. This structure is most 
likely conf used emission from previous activity cycles of 
the BCG dHogan et al. I l2015al ). The integrated emission 
of this amorphous halo component does not account for 
all of the flux in the confused observations. However, such 
high resolution observations at other frequencies to fully 
detangle this amorphous halo emission are not available 
currently. Shown on the SED are the core only observations 
along with highlighted confused, amorphous halo and WAT 
data points. 

RXJ0132-08 

PKS 0130-083, z=0.1489. This source appears to be either 
a low-peaking GPS or a CSS. The source do es not appear 
a bov e the 5<r cut-off of the VLSS catalogue (iCohen et al. I 
120071) but overlaying the VLSS map, the source is present at 
the optical position as a 3cr detection. The central engine is 
resolved in VLBI maps dBourda et al. Il2010l . l201ll : IPetrov I 
1201 ll. [). Fit a GPS model to the integrated emission (that 
returns parameters; 2.55±0.04, -0.24±0.03 and -0.11±0.03) 
and show the VLBI peaks on the SED (although note that 
these are not included in the fit). 

RXJ0352+19 

2MASX J03525901+1940595, z=0.109. The SED is rela¬ 
tively poorly sampled, consisting of only four detections and 
limits from the pointed GIS M O ob servation (at 150 GHz) 
and the VLSS dCohen et al. |[2007i ). A single power-law fit 
returns a=0.41±0.18 and A=20.4±5.8, which is consistent 
with the limits and suggests a flat spectrum source. How¬ 
ever, there could be a small GPS-like component with a tail 
but available data are not sufficient to confirm nor deny 
this, hence the power-law is perhaps the most robust fit 
possible with current data. 

RXJ0439+05 

2MASX J04390223+0520443, z=0.208. The SED can 
be well fit with a strong GPS component, that returns 
parameters 1.42±0.03, 1.67±0.02 and -0.72i0.01 and then 
a steep spectrum tail to lower frequency (less than 1 GHz) 
with a=1.10±0.04 and A=60.2±5.5. 

RXJ0747-19 

PKS 0745-191, z=0.1028. The SED potentially flattens 
below around 200 MHz but other than that presents as 
a steady power-law out to higher than 150 GHz, with 
a=1.26±0.07 and A=3064.7±23.4. 

RXJ1347-11 

GALEX J134730.7-114509, z=0.450. The BCG is sur¬ 
rounded- by a confirmed 500kpc radio mini-halo (see 
iGitti et al. 1120071 ). Using high resolution VLA observations 
from the NRAO Image Archive, an SED of only the BCG 
is created. This suggests an active AGN, whose SED is well 
fit with the GPS model with parameters of Ao=1.35i0.05, 
Ai=0.30±0.01 and A2=-0.35i0.01. The mini-halo emission 
has a steep spectrum with a=1.20±0.04 and A=47.2±4.1. 
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The similarity of the SED of this source to other observed 
spectra in systems without observed mini-halos, but with 
smaller amorphous haloes lends evidence to the belief that 
the seed populations for true mini-halo es are b uil t fro m 
repeated AGN activity in the BCG llHoean et al. 1 12015 al b 


RXJ1350+09 

2MASX J13502209+0940109, z=0.1325. This source ap¬ 
pears to be highly var iable, and has previously been identi¬ 
fied a s a BL-LAC fe.g. iMassaro et al. 1 12009: RichaxdsrfalJ 
l201ll b A non-detection in the VLSS ( Cohen et al. I 2007 1 
shows that this source has no steep component at low 
frequencies and hence the SED appears to be GPS-like. Fit 
the SED with a GPS model, which returns a sharp peaked 
spectrum with parameters Ao=2.42±0.04, Ai=0.70±0.03 
and A 2 =-0.65±0.03. 


RXJ1504-02 

LCRS B150131.5-023636, z=0.2171. This cluster contains 
a confirmed radio mini-halo dGiacintucci et ai\ 1201 l h. 
140kpc in extent. This mini-halo manifests itself as a steep 
spectrum tail in the SED to low frequencies (below about 
1 GHz). There appears to be an active core within the 
central BCG that gives rise to a GPS-like component 
in the integrated SED, that can be fit using the GPS 
model with parameters Ao=1.67±0.01, Ai=-0.09±0.01 and 
A2=-0.16±0.04. The mini-halo e mission can b e iso lated 
from that of the AGN at L-band d White et al. Ifl99~7 ') and 
agrees well with lower frequency data (see SED), giving 
a power-law fit to this lower frequency component with 
a=l.19+0.05 and A=31.4±3.9. 


RXJ1558-14 


PKS 1555-140, z=0.0970. This is an actively evolving 
system with good historical coverage. See section 15.21 for a 
more detailed description of this system. Here we perform 
an average GPS model fit of the core over the previous 
40 years, recovering parameters of 2.69±0.04, 0.31±0.10 
and -0.45±0.06 and then fit a power-law for the steep 
spectrum tail which has a=1.30±0.20 and A=71.6±4.0. 
RXJ1715+57 

NGC6338, z=0.0282. Considering the archival C-band 
measurements (e.g. 1 Gregory fe Condon I 1 199il l suggests 
that this source may show some long-term variability in 
here. Without the GISMO observation the SED would 
perhaps be interpreted as a GPS, however the 150 GHz 
detection shows that the spectrum remains flat up to 
high frequency. The GISMO flux could be interpreted as 
evidence of ‘flickering’ at high frequencies or variability. 
Either way, the SED is most reliably fit using a power-law 
with a=0.41±0.04 and A=53.1+3.1. This shows the source 
to be flat spectrum and is co nsistent with the non-detection 
in VLSS dCohen et al. I [20071 ) . 

RXJ1832+68 

2MASX J18323551+6848059, z=0.205. The SED is domi¬ 
nated by steep spectrum emission to about 10 GHz that is 
fit using a power-law with a=0.92±0.03 and A=194.8±7.7. 
Above this frequency the spectrum appears to ‘bump’ 
outwards. VLBA observations (Hogan et al. 2015c) at 
C-band s how that there is an active core c omponent in here 
(see also lLaurent-Muehleisen et al. lfl997T ). By considering 
measurements of the core flux only, this active component 
can be well fit with a relatively high-peaking GPS com¬ 
ponent with returned parameters of 0.78±0.36, 1.31±0.53 
and -0.56±0.20. When combined with the power-law 


this interpretation appears to well explain the integrated 
spectral shape. 

RXJ2341+00 

PKS 2338+000, z=0.2767. This is an unresolved radio 
source (at few-arcsecond resolution) whose SED is steeper 
at higher frequency yet is flattening at less t han a few G Hz. 
The 365 MHz flux from the TEXAS survey dDouglas et al. I 
1 1996h shows that the SED does turnover. Fit the SED with 
a low peaking GPS, with returned parameters of 2.62+0.01, 
-0.16+0.06 and -0.15+0.04. This source is a no n-detection at 
the 5 a detection limit of the VLSS catalogue dCohen et al. I 
l2007lj , and the limit is consistent with the inferred turnover. 
Looking at the VLSS map overlaid onto an optical image of 
the region, we see that there is a faint (<3<r) detection at 
the correct position. This point was not included in the fit 
but is included on the SED to highlight that it is consistent 
with being a genuine detection. 

Z235 

2MASX J00435213+2424213, z=0.083. The SED of this 
source has a flat spectral index, that can be fit with a 
single power-law of a=0.45±0.06 and A=61.7±7.3. The 
returned index is co nsistent with the non-detection in VLSS 
dCohen et al. IlioOTl b 

Z8193 

B3 1715+425, z=0.1754. This source appears to contain a 
variable co re component that i s conf irmed by VLBI mea¬ 
surements dBourda et al. Il2010l . l201ll . see also Hogan et al. 
2015c, in prep). The SED is fit by a combination of a GPS 
model to account for the spectral flattening above about 
1 GHz and a steep power-law tail to lower frequencies. 
Note that the sourc e is below the dete ction threshold of 
the VLSS catalogue dCohen et al. Il2007fl but appears as a 
«3 <j detection at the correct position in the map overlay. 
The flux is therefore retrived using the AIPS task IMFIT 
from a VLSS cut-out. The returned parameters for the 
GPS component are 1.82+0.01, 0.26+0.01 and -0.26+0.01. 
Accounting for the core contribution, a power law with 
a=0.82±0.04 and A=89.2±5.7 is fitted to the low frequency 
tail. 

Z8276 

2MASX J17441450+3259292, z=0.075. There is a distinct 
split in the spectral index of the SED at approximately 
1 GHz. There also appears to be large amplitude variability 
at higher frequencies. Fit the steeper component using a 
power-law with a=1.26±0.03 and A=80.0+6.0 and then 
fit the GPS model to the core only measurements, which 
returns GPS-model parameters of 1.57+0.50, 0.52+0.76 
and -0.30+0.26. 


This paper has been typeset from a T)5^/ RT(rX file prepared 
by the author. 
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Figure Al. SEDs for sources withing our sample. Typically, data-points represent the spatially-integrated flux density at the given 
frequency. In many cases this is dominated by either the core, or the extended emission, as indicated by the fits. Core only data-points 
have been highlighted for instances where a separate SED was fitted for cores that could be morphologically isolated, to differentiate 
these from instances where the integrated flux shows clear variability. Likewise, other ‘special’ cases have been indicated on the plots. 
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Figure A2. SEDs, continuation of Fig. A1 
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Figure A3. SEDs, continuation of Fig. Al 
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Figure A4. SEDs, continuation of Fig. A1 
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Figure A5. SEDs, continuation of Fig. Al 
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Figure A6. SEDs, continuation of Fig. A1 
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